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CHAPTER I 



INTRODUCTION 

The frequency of discovery of gas condensate reservoirs, 
as well as their importance to the industry, has increased 
greatly in recent years with the trend toward drilling into 
the deeper environment of higher temperatures and pressures. 
Craft and Hawkins ^ point out well test data confirm this 
trend as shown in Figure 1. This graph shows the discovery 
trend for 17 parishes in southwest Louisiana. The reservoirs 
were separated into oil and gas or gas-condensate types on 




Fig. 1. Discovery frequency of 
oil and gas or ga s - c o nd e ns a t e 
reservoirs versus depth. (Craft 
and Hawkins, Applied Petroleum 
Engineering . ) 
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the basis of well test gas oil ratios and the API gravity of 
the produced liquid. While oil discoveries predominate at 
depths less than 8000 f eet ; gas and g as - c o nd e ns a t e discover- 
ies predominate below 10,000 feet. The decline in discover- 
ies below 12,000 feet is due to the fewer number of wells 
drilled below that depth rather than a drop in the occurrence 
of hydrocarbons. 

Many of these high pressure condensate reservoirs are 
discovered in conditions above their dew points. With pres- 
sure reduction, when placed on production, they allow retro- 
grade condensation of the valuable heavier hydrocarbons within 
the reservoir rock. Once condensation has taken place within 
the reservoir, this liquid cannot be produced without con- 
siderable expense. The mechanisms by which this condensation 
may be retrieved is either by r e v apo r i za t io n or by repressuri- 
zation by injected dry gas. However, several pore volumes of 
dry gas must usually be injected before sufficient revaporiza- 
tion occurs, making this process costly. Even then, the dry 
gas may not sweep all the volume in which condensation has 
taken place. In order to retrieve this condensate, which is 
many times more valuable than the gas itself, or to prevent 
its condensation within the reservoir, it is necessary to in- 
vestigate the possibility of dry gas cycling and pressure 
ma i n t e na nc e . 

Cycling of the dry residue gas may act to maintain the 
pressure in the reservoir at or near the dew point, thereby 
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preventing the occurrence of the retrograde condensation 
phenomena. It may also tend to revaporize the condensate and 
sweep it, along with the cycled gas, out of the reservoir. 

Even though it is apparent that cycling a reservoir 
which is susceptible to the retrograde phenomena will produce 
more of the valuable hydrocarbons in place, the economic fac- 
tors involved are all important in deciding whether to invest 
in this procedure. Roughly, a reservoir that contains as much 
as fifty barrels of condensate per million standard cubic feet 
of gas is a good prospect for gas cycling. On the one hand, 
the dry gas sales will be postponed until the cycling is prac- 
tically complete, resulting in decreased present value of the 
gas. However, it is generally possible to produce the con- 
densate faster, hence earlier, when cycling is adopted and 
more dollars are generally obtained due to the increased ulti- 
mate recovery. The purpose of this thesis is to evaluate whether 
or not these advantages are worth the greater investment in- 
volved for the field studied. 

Data are obtained from an actual field on the Gulf Coast 
of Louisiana containing approximately 600 billion SCF of gas. 

The size of the reservoir was determined from geologic data 
and from production data to date. Future prediction of reser- 
voir performance was made for pressure depletion, normal cycling 
and pressure maintenance methods of exploitation. An economic 
comparison of the cases studied is presented. 



CHAPTER II 



DISCUSSION 

A . Nature of Gas Condensate Reservoirs 

Gas reservoirs are defined as having a lean gas, con- 
taining a minimum amount of the heavier hydrocarbons, where 
the composition of this produced gas does not change with 
pressure reduction as the reservoir is depleted. Oil reser- 
voirs range in gas content from zero, or dead oil, to a gas- 
oil ratio of a few thousand cubic feet per barrel. In gas 
reservoirs, on the other hand, the liquid may be vaporized 
into 100,000 cu. ft. per barrel or more. Gas condensate 
reservoirs are classified between these two extremes. Gas 
condensate reservoirs, or distillate reservoirs as they are 
sometimes called, usually produce light colored or colorless 
stock tank liquids with gravities above 45 degrees A.P.I. at 
gas-oil ratios in the range of 3,000-100,000 SCF/bbl. The 
production from these reservoirs is predominately gas from 
which some liquid is condensed in the surface separators, 
which phenomenon gives them the name gas - condens a t e . 

Another definition of the three types of fields described 
can be derived from an interpretation of the p r e s s u r e- t emp e r a tu r e 
phase diagrams. (See Figure 2 for a typical diagram.) When- 
ever a reservoir fluid is discovered at a temperature above the 
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c r i c o n d e n t h e r m , or maximum two phase temperature for the fluid, 
the composition of the produced fluid will not change as the 
reservoir is depleted since the reservoir temperature can be 
expected to remain fairly constant (Path A-A^ on the diagram). 




Fig. 2. Pr es s u r e - t emp er a t u r e phase 
diagram of a reservoir fluid. (Craft 
and Hawkins, Applied Petroleum Reservoir 
Engineer ing . ) 



There may, however, be liquid produced as the fluid flows up 
the well bore and into the separators if temperature reduction 
below the c ircondentherm occurs within the tubing string (Path 
A-A on Figure 2). If, on the other hand, the reservoir is 
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found to be at a temperature below the critical temperature 
(T^) and pressure above the bubble point, the reservoir would 
be in a one phase liquid state. This is called a bubble 
point reservoir and the fluid will remain a liquid although 
as the bubble point pressure is reached, with pressure deple- 
tion in the reservoir, gas will come out of solution in the 
oil and will be produced along with the oil in an ever increas 
ing quantities (Path C - C ^ in Figure 2). In dew point or retro 
grade gas condensate reservoirs the reservoir fluid exists at 
a temperature between the critical temperature and the cricon- 
dentherm and at a pressure above the dew point pressure of the 
fluid. In this situation, when the pressure declines, composi 
tion of the reservoir fluid, which may be termed gas in this 
state, will remain the same until the dew point is reached. 

At this point as the fluid expands i s o t h e rma 1 1 y , a liquid be- 
gins to condense out of the gas and will adhere to the walls 
of the rock, tubing, etc., leaving a gas leaner in the heavier 
hydrocarbons to be produced. The term retrograde is given to 
this phenomenon since one would expect vaporization to take 
place with a pressure decline rather than condensation. At 
some intermediate pressure, between the dew point and abandon- 
ment, this retrograde phenomenon will reach a maximum, below 
which pressure the condensate will revaporize into the gas. 

( Pr es s ur e- Temp era tur e phase relationships for g a s - c o n d e ns a t e 
reservoirs are shown on Figure 2 by the path from B to . ) 
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B . Retrograde Loss of Condensate 

Liquid that is condensed inside the reservoir rock 
remains immobile in the usual case and therefore is lost to 
production unless or until revaporization can take place 
either by further pressure decline, repressurization, or by 
sweeping with a dry gas. The alternative to revaporization 
of this condensed liquid is to prevent it from condensing in 
the first place and this possibility leads to the investiga- 
tion of pressure maintnenace by injection of dry gas into the 
reservoir. 

C . Experience in Gas Cycling 

In general, the gas with the greater condensate content 
exhibits the greater condensate loss upon pressure depletion 
so usually it proves to be the more profitable reservoir to 
cycle with gas. In gas cycling the condensed liquid is sepa- 
rated from the wet (produced) gas in a gasoline plant or other 
surface facilities. The residue or dry gas is then re-injected 
into the reservoir where it maintains the reservoir pressure 
and also serves to drive the wet gas toward the producing wells. 
This procedure is continued until the dilution by dry gas at 
the producing well renders the operation unprofitable at which 
time the reservoir is blown down, or allowed to produce nor- 
mally through pressure depletion. During this pressure blow- 
down not only will the lean gas in the reservoir be produced 
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but also additional liquid may be recovered from both the 
dry gas invaded zone and the uninvaded portion of the 
res ervo ir . 

2 

Bowers has shown that cycling in the Grapeland Field, 
Houston County, Texas, resulted in a recovery of seventy-five 
per cent of the initial butanes and heavier hydrocarbons in 
the reservoir, while straight pressure depletion would have 
recovered only thirty-nine per cent of the initial butanes-- 
plus. The liquid that could have been recovered by pressure 
depletion, then, was only fifty-two per cent of that which 
was actually produced by the cycling process. In this example 
it was shown that the ultimate liquid recovery was reached in 
approximately eleven years but twenty-one years would have been 
required by straight depletion, thereby making the present worth 
of the liquid products much greater with gas cycling. 

In a study of the Krotz Springs Field, St. Landry Parish, 

3 

Louisiana, Robertson predicted that twenty per cent more con- 
densate could be recovered by cycling, resulting in a gross 
cash gain after taxes of $2,000,000 for the twenty-five year 
period studied. 

4 

T. W. Brinkley concluded that benefits in improved 
recovery of stock tank condensate by pressure maintenance 
methods may be as great as 300 per cent in rich condensate 
reservoirs but grade progressively downward to negligible 
benefits for the lean condensate reservoirs with GOR' s 
greater than 300,000 cubic foot/bbl. 
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Although cycling may appear to be the ideal solution 
to the retrograde condensation problem there are other fac- 
tors that may make this procedure uneconomical. The deferred 
income from the dry gas sales may prove to be a significant 
factor in the economic analysis of the field. Another eco- 
nomic consideration is the added investment required for 
additional injection wells, gas compression and distribution 
system to return residue gas to the wells and for a liquid 
recovery plant. Also, it must be remembered that even with 
gas injection at a pressure above the dew point all of the 
liquid hydrocarbons may not be recovered. 

Three different recovery factors must be applied with 
gas cycling. The microscopic displacement efficiency is 
approximately 70 to 90%. The volumetric sweep efficiency, 
or the per cent of initial pore volume invaded by the sweep- 
ing gas to abandonment of the producing wells can vary from 
50 to 907 o . Finally a permeability stratification factor 
must be applied to take care of the problem of gas sweeping 
through the more permeable strata and reaching the producing 
wells before the tighter stringers can be swept. The over- 
all condensate recovery factor for a gas cycling project is 
the product of these three factors. On the plus side of 
recovery by cycling is the fact that additional liquid is 
produced during blow-down from the less permeable or unswept 
portions of the reservoir as these products bleed into the 
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more permeable streaks and ultimately to the wells. Also the 
efficiency of extracting condensate from the wet gas is higher 
with the type of separator plant that is generally used in a 
cycling project. 

Any reservoir exhibiting retrograde properties and hav- 
ing a richness of 50 bbls/MMSCF of condensate is a good pros- 
pect for cycling of the dry residue gas to improve recovery. 
Each reservoir, however, is an individual case and, as such, 
requires a separate study. Variables such as permeability 
stratification, reservoir geometry, gas properties, and market 
availability can team up to make one reservoir an extremely 
profitable enterprise while another may prove entirely infea- 
sible. This thesis presents the study of one field based on 
the fluid properties and geologic structure of that one field. 



CHAPTER III 



DEFINITION OF THE FIELD STUDIED 

A . Development of the Field 

Production from this field began in January, 1961 and 
from that time two to three wells have produced continuously 
as shown in Table II. The pool is one in which higher than 
normal pressures are encountered. It has been obvious from 
the low gas-oil ratios measured that the produced gas was 
rich in condensate. By January 1, 1964, a total of 20 billion 

standard cubic feet of wet gas had been produced from the two 
wells completed in the zone studied. The pressure had dropped 
from an estimated initial pressure of 8838 psia to 8100 psia 
as shown in Figure 7. Based on these production data, calcu- 
lation of the initial gas in place yielded the value of 606.9 
billion standard cubic feet. These values form the basis for 
future performance predictions from the base date of January, 
1964. 

B . Structure and Reserves 

Interpretation of available geologic data indicates a 
deep-seated salt domal structure. Producing depths at which 
this study was conducted were from approximately 11,300 to 
11,600 feet. A structure contour map of the top of the pro- 
ducing formation is presented as Figure 3. The productive 
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interval is taken to be a shaly-sand with intergranular 
porosity. 

There are a total of nine wells drilled through the 
formation studied and from which data are available from 
various kinds of logs. Porosity was determined using micro- 
logs from some wells and sonic logs from others. The porosity 
throughout the formation is fairly constant, varying from 25 % 
to 297, in the pay zone. Although shaly streaks were noted, 
this fact was taken into consideration in counting sand, 
wherein the shale sections were discounted. For purposes of 
the material balance calculations an average uniform porosity 
of 2 87, was used. 

Electric logs or indu c t io n- e 1 ec t r i c logs were available 
from all nine wells. Three of these logs are shown in Figure 4 
These logs were used to compute formation water resistivity 
using the Spontaneous Potential Curve. Water saturation was 
calculated from the induction log or lateral curve as avail- 
able. Water resistivity was found to be approximately .03 ohm- 
meters at formation temperature and this value w,as used through 
out the field. Computed water saturation tended to vary from 
8-1 07o near the top of the gas zone to about 2 57, nearer the 
bottom of the pay and it grades to 10 07, water saturation in 
those wells where a gas-water contact is observed. A value of 
2 07o wa t er- s a tur a t ion was used for the entire formation as being 



the weighted average value. 



WELL LOG CORRELATION 
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Figure 
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Side wall cores were available from one well. The 
porosity values obtained from the sidewall cores checked pre- 
cisely the values obtained from the sonic logs and micrologs. 
The permeability, however, did not appear uniform, nor did 
these values check well against the effective permeability 
calculated from a pressure draw-down test or against permea- 
bility calculated from a resistivity gradient observed in some 
of the wells. The six values obtained from side-wall cores 
were 84, 134, 50, 512, 131 and 116 mi 1 1 i da r c i e s . The effec- 

tive permeability to gas obtained in the pressure draw-down 
test proved to be 50 millidarcies and this was the value used 
in the material balance calculations since the method of tak- 
ing side-wall cores tends to result in fractured, non- 
representative samples of the reservoir. A plot of the data 
taken from the pressure draw-down test is shown in Figure 5. 

For volumetric sweep-out calculations, 50 millidarcies 
was used for permeability while the effective thickness of pay 
was allowed to vary over the field. This gave a value of flow 
capacity for each point in the reservoir which was used in the 
sweep efficiency program. 

Connate water saturation was found to be 20 per cent. 
This was an average value based on water saturation in the 
lower portion of the gas zone. Residual gas saturation to 
water is estimated to be 30 per cent. This is justified by 
the gas saturation near the top of the water column which, it 
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Pressure Drawdown Test 



Figure 5 
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can be assumed, has encroached into the reservoir to some 
extent during a recent producing time. 

Determination of net effective thickness in this forma- 
tion was very difficult because of the shalyness of the sand 
and because of the meager geologic information. Only nine 
wells have penetrated the 2000 acres of estimated pay forma- 
tion. The locations of the two faults--one thought to have a 
closure of about 100 feet in the northern extremity of the 
pool, the other with an estimated closure of about 120 feet 
in the south--are by no means exactly defined. An error of 
the position of the fault to the north of only a few feet 
would produce a substantial error in the estimated reserves 
because of the great thickness of pay adjacent to this fault. 
Calculation of reserves based on the volumetric method is 
presented in. Table II. The east and west boundaries of the 
field are located at a gas-water contact which has shown 
little or no encroachment during the current gas production. 

An isopach map produced from the nine control wells is 
presented as Figure 6. 

A volumetric gas balance was made using the production 
data to date and bottom hole pressure tests taken in 1961 and 
1963. A bottom hole pressure was obtained from the draw-down 
and build-up pressure test conducted in connection with this 
study. This gas balance was calculated on a Control Data 
Corporation 1604 high speed computer using a program developed 
by the Union Oil Company of California. 
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J. Donald Clark, Chief Reservoir Engineer for the Gulf 

Division of the Union Oil Company of California, has developed 

and applied or supervised the application of a very simple 

transient energy balance equation to numerous natural water 

12 

drive gas reservoirs found in the Gulf Coast area. A primary 
requirement of this equation is the derivation of the gas for- 
mation volume factor, B , to convert standard cubic feet of 

g 

gas to reservoir barrels. This conversion for any particular 
reservoir is presented mathematically as follows: 



15.025 Tf x _Z 0051505 Tf 2/p 
520(5.61) 



Where: 15.025 is base pressure for Louisiana (psi) 

520 is standard temperature °R 
5.61 converts cubic feet to barrels 
T is reservoir temperature, degrees Rankine 
Z is gas deviation or compressibility factor 
P is the reservoir pressure, psia 
i as a subscript will denote initial or 
original conditions 



The cumulative reservoir barrels of gas produced or withdrawn 
is equal to the cumulative reservoir barrels of gas expansion 
in the reservoir, or 



2 . G B 
pw g 



G ( B -B 
g 




Where: Gp W is cumulative standard cubic feet of 

wet gas produced 

G is standard cubic feet of original gas 
in place 



In the case of water drive reservoirs and water production 



the statement of the equation then becomes: (Cumulative Net 



Reservoir Withdrawals) - (Cumulative Reservoir Expansion) 
(Cumulative Water Influx), or 
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3. (G B + W ) - G (B - B . ) = W 
pw g p g gi e 



4. (G B + W ) - G ( B - B . ) = AEAPQ 
pw g p g gi r 



The water influx equation, W^ = AEAPQ^, is a simplified 
method of writing the equation presented by Hurst^ in his 
paper entitled "Water Influx into a Reservoir and Its Appli- 
cation to the Equation of Volumetric Balance." Details of 
solving the water influx equation were presented by van 
Everdingen and Hurst. The "A" is the water influx constant 
which encompasses conversion factors, porosity, radius of 
reservoir, effective reservoir thickness, effective compressi- 
bility of reservoir fluid, viscosity of reservoir fluid and 
fraction of reservoir perimeter exposed to water drive. The 
above parameters are constant for any particular reservoir 
and when used with the volumetric balance equation can be 
solved as an all encompassing constant for that particular 
reservoir. 

The true simplification is in the manner of arriving at 
the value of EAPQ^ for various balancing periods. 

2 

The dimensionless time equation, t^ = KT/^C^r^ is 

2 

derived from the hydraulic diffusivity factor, a = K/p^C 

13 

presented by Hurst in 1933. The K, or permeability factor, 
actually is the one most difficult to measure since it is 
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effective permeability to water, which we seldom if ever 
obtain. This value could be found by running draw-down tests 
in salt water wells, but actually it can be estimated to a 
sufficient degree of accuracy. Viscosity of water can be 
obtained from various sources to a reasonable order of magni- 
tude. Therefore, it becomes important only to establish an 
order of magnitude value of t^, in order to extract reasonable 
values of the corresponding Q , the dimensionless water influx, 
from the tables. 

Experience has shown that a proper time period should 
not exceed three months or approximately a 90-day producing 
period when attempting to strike a balance, especially if one 
plans to calculate future predictions. The constant "A" be- 
comes an over-all correcting constant of transient flow, and 
allows us to solve very simply a satisfactory value of original 
hydrocarbons in-place and the total water influx at any time in 
the producing life of the reservoir. 

Cumulative Production and Pressure History 

Production from the reservoir is tabulated on a monthly 
rate and cumulated on a three-month basis. These data are 
broken down between oil, gas, condensate and water (both 
fresh and salt water). When rates of production have consider- 
able variation, a monthly production rate should be plotted as 
shown on Figure 7. This plot is important in order to show 
the effect of rate on pressure when natural water drives exist. 
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Gas Produced MMMSCF 



Production History 
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Bottom hole pressure measurements are plotted versus 
time on the same graph as production rate. The length of 
shut-in time prior to each pressure measurement should be 
approximately equal and the individual pressures should have 
reasonably close values. If not ; build-up or drawdown tests 
should be run to determine the cause of pressure variation. 

A pressure curve is then constructed through these evaluated 
pressure points. Experience and ingenuity are of tremendous 
value in constructing a reasonable curve. The production 
rate graph must be studied in order to prepare the most logi- 
cal pressure- time curve. Pressure for uniform periods of pro- 
duction are taken directly from this graph. 

Water Influx Data and Basic Calculations 

The next step requires the gathering of basic informa- 
tion to approximate dimensionless time data and corresponding 
dimensionless rates. 



KT (86,400) 

u0 ( 0 . 5 X 10 -4 ) r b 2 (930) 



Here geological data can be used to help determine the radius 
of the hydrocarbon reservoir. Only a close approximation is 
necessary. Compressibility of salt water can be found in the 
literature. Viscosity of salt water can be approximated by 
using chemical handbooks ; available texts, and various publi- 
cations; porosity is approximated from core analyses or logs 
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and permeability to water can be estimated from air permea- 
bility of cores analyzed. 

Dimensionless time values are calculated for cumulative 
three-month periods. The corresponding values are taken 
from the table. Integration of the values for AP times 
over the cumulative time intervals yields the value which, 
when multiplied by the constant A, provides the cumulative 
wa ter influx. 

Gas Conversion to Reservoir Barrels 

The vapor equivalent of the distillate and fresh water 
produced must be added to the amount of dry gas produced. The 
constant for converting liquid water to the vapor state is 
7368.48 cubic feet per barrel. The distillate conversion fac- 
tor was calculated from the hydrocarbon analysis of the sepa- 
rator liquid shown on page 89, Appendix II. This factor was 
found to be 768 standard cubic feet of vapor per barrel of 
distillate. 

Method of Least Mean Squares for Striking an 
Energy Balance 

The method of least mean squares is used for striking 
an energy balance for the natural water drive gas reservoir. 

(G + W ) - G ( B - B . ) = A (EAPQ ) 
pw ps' g gr r 7 

The two unknowns are G, the original wet gas in-place, 
and A, the water influx constant. The method herein presented 
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requires the solution for single values of A and G to meet all 
periods of past producing time. We, therefore, can arrive at 
two equations with two unknowns and place proportional weight 
to the length of production time: 

[Z(G pw + W ps )(ZAPQ r )] - G[(B g - B g .)(ZAPQ r )] = A[(ZAPQ r ) 2 ] 

[ Z ( G + W ) ( B -B.)]-G[(B -B.) 2 ]=A[( ZAPQ ) ( B - B . ) ] 

pw ps' g gi /J g gi' J x r' v g gi' J 

These two equations can then be solved for A, the water 
influx constant, and G, the original wet gas in-place. 

In order to check the validity of the values of A and 
G obtained by this method, multiply the LMSA value of A times 
ZAPQ r and solve for the apparent value of G. This is called 
"periodic check of balance for G assuming LMSA values of A 
as constant," The van Everdingen and Hurst solution to the 
diffusivity equation is here applied to single-phase fluid 
transient flow into a well bore for determination of the rate 
and volume of encroachment of water. The gas field is assumed 
to be embedded into an aquifer of large extent. In this 
instance, the reservoir is taken to be a large well, the 
radius of which is the equivalent field radius, and the aqui- 
fer is considered to be the reservoir from which fluid is 

g 

withdrawn and water is expelled. 

The computer program, indicating method of computation, 
reservoir conditions, initial gas in place, and cumulative 
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water influx from past production data, is included as 
App end ix I . 

Since the available geologic information is considered 
isufficient to describe the reservoir accurately enough for 
calculation of reserves, the results obtained from the volu- 
metric gas balance were used for estimating original gas in 
place and water encroachment. These results were also the 
basis for the future performance prediction and economic 
analysis of alternate methods of exploitation. The original 
gas in place was found to be 606.9 billion cubic feet measured 
at standard conditions (15.025 psi at 60°F). Water influx was 
found to be negligible. 

C . Reservoir Fluid Characteristics 

Samples of both the separator liquid and gas were ob- 
tained in March, 1964, and recombined in accordance with the 

measured gas-liquid ratio (stock-tank) of 17,150 SCF separator 

* 

gas per bbl. stock tank liquid. It was difficult to obtain 
constant gas-oil ratios for the purpose of getting a repre- 
sentative sample and flowing at a constant rate for about three 
days was required before a steady GOR was achieved. Analysis 
of the recombined sample indicates a dew point pressure above 
the present bottom hole pressure. Retrograde condensation 
was taking place and two phase flow occurring near the well 
bore caused the fluctuating producing gas-liquid ratio. Reser- 
voir temperature at time of sampling was 240 degrees farenheit. 
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Primary separation temperature was 125 degrees farenheit and 
primary separator pressure was 1020 psig. Dew point pressure 
was found to be 8260 psia and gas richness was 66.07 bbls. of 
stock tank liquid per million standard cubic feet of well 
stream effluent. Detailed results of the P-V-T investigation 
are included as Appendix II. 

Samples of separator gas and separator liquid were ana- 
lyzed by the conventional low temperature distillation method. 
The results are reported showing both the composition of each 
sample and the computed analysis of the well stream based on 
the gas-liquid ratio of the primary separator. The separator 
liquid production was calculated from the measured stock tank 
production by applying the determined shrinkage factor. 

Following the compositional analyses, portions of the 
primary separator liquid and gas were physically recombined 
in their produced ratio in a variable volume, g 1 a s s - w i n dow e d 
equilibrium cell. Determinations on this mixture were divided 
into the following two main categories: 

I. Dew point pressure determination and pressure- 
volume relations on a constant weight of reservoir fluid 
at the reservoir temperature: The procedure consisted of 

establishing equilibrium between gas and liquid phases at 
a low pressure and measuring the volumes of liquid and 
gas in equilibrium at that pressure. The pressure was 
then raised by the injection of mercury into the cell and 
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phase equilibrium established again at the higher pressure. 
This procedure was repeated until all of the liquid phase 
had vaporized, at which point the saturation pressure was 
observed. The cell pressure was then raised above the dew 
point pressure in order to determine the sup er compr es s i- 
bility characteristics of the single phase vapor. As a 
check on all readings and particularly to verify the dew 
point, the cell pressure was incrementally reduced, equi- 
librium established and volumetric readings made. Reported 
in Appendix II are the relative volume relations and spe- 
cific volumes of the reservoir fluid over a wide range of 
pressures as well as deviation factors of the single phase 
vapor above the dew point, (Also reported in Appendix II 
are the dew point pressures resulting from recombinations 
at gas-liquid ratios above and below the ratio measured 
at the time of sampling.) 

II. Compositions of the produced well stream and the 
amount of retrograde condensation resulting from a stepwise 
differential depletion: This procedure consisted of a 

series of constant composition expansions and constant vol- 
ume displacements with each displacement being terminated 
at the original cell volume. The gas removed during the 
constant pressure displacement was charged to the low 
temperature fractional distillation equipment for analysis. 
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D . Conf ormanc e 

The description of the simultaneous flow of fluids 

through porous media in terms of relative permeability and 

capillary pressure has been adequately discussed in the lit- 

9 

erature. For example, Muskat in Chapter VII briefly dis- 
cussed the dynamical effects associated with capillary 
phenomena. He also outlined the theory of po t e n t iome t r i c 
models and illustrated the application of the po t en t iome tr i c 
model to cycling systems. Douglas, Peaceman, and Rachford,^ 
presented a method for calculating mu 1 1 i - d ime ns io na 1 immiscible 
displacement based on a numerical solution of a finite differ- 
ence analogue describing the displacement process. This tech- 
nique insures that at any point in the sand, the flow of oil 
per unit area across the direction of flow can be represented 
by a vector which is given by 

K 

Tf = -K — (Vp + p g Vh). 
o u o o 



The operator, V, equals ~t + T ; 



j are unit 



)x oy 

vectors in the x and y directions, respectively; k is the 



local permeability, k^ o is the relative permeability to oil, 

(i Q the oil viscosity, p^ the pressure in the oil phase, p Q the 
oil density, and g the acceleration of gravity. Similarly, 



k 

u = -k — — (VP + P g vh) 

W (! W W 



Fluid densities are assumed constant and oil and water 
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potentials are defined 




<5 = p + p gh 

x w w w 



The author proceeds to solve the problem by a system of dif- 
ference equations. 



Dr. S. J. Pirson at The University of Texas is followed. In 
this method the reservoir is divided into a grid and each 
intersection of the grid is treated as a point (I, J) within 
the reservoir where Darcy's law may be applied. The flow 
across the grid between these points ; q, is now analyzed. 

The rate of flow over a distince Ax across an area open to 
f low, Ay, can be determined by Darcy's law: 



In solving this reservoir a technique demonstrated by 








J 





q = 1.12 7 Ax [ 

1 



P 




o 



Ay 



In order to satisfy the continuity equation the vector 



sum of the flow in all four directions must equal zero unless 



there is fluid withdrawal or fluid injection at the point 



Q = q L + q 2 + q 3 + q 4 



A = 



1.127 



[K. (P . . _ * . . 

h x i-l. J l.j 
av ; J 



P, . )+K, (P, ,.,-P, . ) + K, (P.^. ,-P. 



h ' ‘ i, j + 1 ‘ i, j 

a v 7 ' 



a v 



i+l, j 



+ K, (P. . ,-P. .)] 

h i , J - 1 i , j 
av ’ J ’ J 



It is seen that without injection or production from the 
point studied Q=0 and the Kh is the average flow potential 



between adjacent points. This equation is now solved for ^ 

and we find: 

2Q3n„ 



+[ K , .+K . , )P 



+ (K, . . +K . )(P, , , ) + 



1.127 i+l,j i,j i+l,j i-l,j i , j i-l,j 



P . 

1 > J 



+(K. . +K . .)(P. . - )+(K . . +K . ,)(P. . .) 

i,j+l i,j /v i,J+l i,J-l i,J i,J-l 

K. ,,, + K. . + K, , , . + K. . . + 4K. . 

i,J + l ijJ-1 i+l^J i - 1 , J 



The boundary conditions are those restrictions on the 
limits of the connected region and the position and magnitude 
of the imposed flows or pressures. At the boundary of the 
permeable sand, the velocity of flow perpendicular to the 
boundary is zero, so it is sufficient to establish equal poten 
tials across the boundary. Since (V(Jj) = 0 across the boundary 
no flow occurs. This same condition occurs at a fault bound- 
ary. As the water saturation increases near the water table, 
the relative permeability to water increases greatly while 
that to oil diminishes. By assigning these relatively high 
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values of flow potential in the equation aquifer pressure is 
maintained in the advancing flood front and movement of the 
front is regulated across this boundary. 

In solving this potential equation for an oil reservoir, 
the term -2(3^/1.127 is treated as a constant. However, 
since the flow equation depends on pressure, the constant for 
a gas field is augmented by a variable pressure term: 



q =-1.127 



K dP 
\i dx 



(q)( p )(T)( z ) 

sc sc 

T5 . 6T5)TT s " c )(A) (¥) 



For Oil 



For Gas 



Volumetric sweep efficiency was investigated by this 
technique with the aid of a mathematical model which was run 
on the CDC 1604 high speed computer. This : computer program 
is included as Appendix III. The program assumes steady-state 
flow at a point, and sequentially at each point in the reser- 
voir in order that Darcy' s flow equation may be applied. Using 
an alternating direction iteration procedure, the computer then 
reduces the total flow at each point in four directions to zero 
and by many iterations establishes each potential from the 
adjacent four potentials in order to achieve a balance. At 
each point selected on the model the quantity of flow into or 
out of the reservoir is assigned in the event a well is located 
at that point. The flow potential, or the product of effective 
permeability times the thickness of effective pay, is also 
assigned to each point. By applying Darcy’s law in each of 



33 



four directions, by summation of these directional flow rates 
to zero, and by solving this equation for the potential at 
each point, a potential distribution map of the reservoir is 
produced. Isopotential contour lines are now drawn on a scale 
map of the reservoir. The primary criterion for equivalence 
between the model and the actual reservoir is the creation of 
geometrically similar potential fields which are determined 
only by the thickness, permeability, and boundary conditions. 
Flow lines were superimposed on the potential map and the 
method of squares utilized by Seba^ is used to define the 
flow tubes through which a constant quantity of reservoir 
fluid would be expected to flow. 

Flow lines of a fluid in a porous medium, as defined 
in Darcy's law, are everywhere orthogonal to the isopotential 
lines. The method of squares employs the fact that if a family 
of curves is constructed such that they are orthogonal to one 
member of another family of curves, then the second family of 
curves may be drawn orthogonal to the first. This is done by 
constructing them such that adjacent distances between the 
streamlines are equal to adjacent distances between the iso- 
potentials. A stream tube is a two-dimensional flow system 
defined as the area between two consecutive stream lines. Flow 
rate can be defined by the following: 

Q = <h)(- ^XAxH-ff) 



a v 
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In any one square Ax = As and Q is constant so: 

2 

At = (constant) 

The increment of time can now be calculated for each of the 
squares within each flow tube so that the position of the 
front in each flow tube can be drawn at any selected time 
period. This procedure is repeated for various injection- 
production well combinations until the maximum volume of the 
reservoir is swept. Rearrangement of the well positions in 
the reservoir matrix is achieved simply by changing the quan- 
tity of flow into or out of that point in the reservoir. Sweep- 
out patterns were predicted with one, two and three injection 
wells and by using four producing wells. The locations of 
existing wells were used as much as possible but alternate 
new locations were also tried. The optimum combination was 
found to be two injection wells and four producing wells. 

Wells marked "Well 4" and M Well 10" on Figure 8 are to 
be recompleted as injection wells. Two new wells are required 
at the locations marked ,f Well 5” and "Well 11. " An existing 
well at location "Well 9 " is to be recompleted in the pay zone 
studied since it is now producing from a deeper zone. "Well 6" 
is producing from this formation and will continue to be used 
for this project. 

A machine print-out of the potentials has been reduced 
to a map shown in Figure 8 picturing the isopotential contour 
lines and on which the frontal position has been drawn at one 
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time just after breakthrough of the injected residue gas. The 
volume of the reservoir swept at this breakthrough time was 
found to be 65%. 

The ability of the dry absorber gas to vaporize the 
condensate in the sand was determined by pressure-volume- 
temperature tests run on a recombined sample of the gas. 

Results of these tests are given in Appendix II. The overall 
displacement efficiency of the residue gas was found to be 
92%. 

The effect of bypassing on a local scale, due to permea- 
bility variations, was analyzed as demonstrated by Standing, ^ 
by arriving at a permeability variation for the sand. Although 
insufficient data are available to establish firmly the perme- 
ability variations throughout the reservoir, a variation .30 
was assigned. This value gives a sweep efficiency at break- 
through of 87%. 

Using the product of these three conformance factors, 
the overall volume of the reservoir swept out at breakthrough 
was 51%, with the dry gas cut after breakthrough increasing in 
the manner calculated by Standing^ in the paper mentioned before. 
Cycling of gas was discontinued when proportion of reservoir gas 
produced was reduced to 61 % by dilution with cycled residue gas, 
due to the cycling economic limit having been reached. 



CHAPTER IV 



PERFORMANCE PREDICTION 

A . Pressure Depletion Method 

Past production has averaged 25 million standard cubic 
feet per day from the reservoir studied. Until this time no 
firm estimate had been made on the reserves nor had any analy- 
sis been made to predict the future performance and probable 
life of the reservoir. From the past production figures the 
initial gas in place was estimated and the pressure decline 
and water influx was analyzed. From these data it is possible 
to predict the future performance of the field. This was done 
through the use of the computer program included in Appendix IV. 

Because of insufficient pressure decline and the lack of 
back pressure data it is not possible to predict d e 1 i v er ab i 1 i t y 
deterioration with declining pressure. It is anticipated that 
as liquid saturation around the well bore increases with static 
pressure reduction and with more severe pressure drawdown, 
deterioration of de 1 iv er ab i 1 i t y will result. For the purposes 
of this study, in absence of definitive data, the rate of pro- 
duction of dry gas was assumed to remain constant and at the 
current rate. The water influx constant obtained from the 
volumetric gas balance program was used to obtain future per- 
formance prediction as well as the calculated quantity of 
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initial gas in place and the adjusted pressure decline 
f igur es . 

The amount of condensate produced, by year, was calcu- 
lated from the P-V-T data obtained in the material balance 
calculations performed on the recombined gas sample placed 
in the equilibrium test cell. In this cell the pressure was 
depleted by removing successively small quantities of gas from 
the cell. Thermodynamic equilibrium was established, material 
balance calculations were made, and retrograde liquid volumes 
and the gravity of the effluent gas were measured at each 
step. The gas removed in each depletion step was charged to 
the low temperature distillation equipment for analysis. The 
produced volumes of reservoir gas, the determined compositions, 
computed gallons per thousand standard cubic feet (G PM) con- 
tent and respective deviation factors are presented in Appen- 
dix II. The volume of retrograde liquid resulting from this 
depletion study is also shown in Appendix II, both in terms 
of barrels of reservoir liquid and per cent of hydrocarbon 
pore space. Shown also are the compositions of the gas and 
liquid remaining in the reservoir after depletion to abandon- 
ment pressure. The amount of condensate calculated to be pro- 
duced with the well stream at each pressure decrement in the 
future prediction was tabulated. The value of gas sold under 
the terms of the present contract is $.20 per MCF and conden- 
sate is sold for $2.95 per barrel. These prices are assumed 
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to be firm for future sales. Operating costs are averaging 
$.00423 per MC F produced and product transportation charges 
for the condensate sold are averaging $.1375 per barrel. 
Application of these values to depletion of the reservoir 
over a fifty-year period with net income discounted at 6 % 
results in a present worth of $38,749,249. It is seen that 
in this case of pressure depletion to 1500 psi only 73.51 per 
cent of the initial gas in place was produced. Since the gas 
is to be delivered to a pipeline at a pressure in excess of 
1000 psi, further recovery would require recompression to the 
sales gas line. This secondary recovery process was not con- 
sidered in this study because the value of gas recovered 50-odd 
years hence is not relevant for a study of possible investments 
now . 

It is significant to note that the condensate recovered 
by pressure depletion amounts to 15.4 million barrels which 
is scarcely 38% of the more than 40 million barrels of recov- 
erable oil initially in place. 

B , Pressure Maintenance through Gas Cycling 

A program of gas cycling was considered wherein gas 
sales from the reservoir would be terminated and all gas pro- 
duced from the reservoir, together with enough make-up gas to 
allow for compression and other gas losses, would be reinjected 
after extraction of all condensable heavier hydrocarbons. As 
determined by the sweep efficiency model studies, two wells 
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should be recompleted as injection wells and four wells would 
be used as producing wells as shown in Figure 8. Since the 
installation is located off shore the wet gas would be pro- 
cessed in low temperature units on platforms near the produc- 
ing wells. The residue gas would be piped from the separators 
to the injection platform which would be located at one of the 
injection wells. Double extra heavy pipe would be necessary 
to carry the compressed injection gas to the other injection 
wells since injection pressures in excess of 8000 psi will be 
encoun t er ed . 

Condensate produced in the separator units would be 
metered and then reinjected into the gas sales line of the 
gas company where it would be transported with other wet gas 
to the on-shore extraction plant for processing. 

Under pressure maintenance it is assumed that the reser- 
voir will continue to produce a gas of the richness measured 
in the most recent flow tests until breakthrough of the in- 
jected dry gas into the first producing well. From this point 
in time the produced gas would decline in richness in direct 
proportion to the residue gas cut produced. When the amount 
of injected dry gas obtained through the production wells 
reached 1/3 of the total gas processed it was found uneconomical 
to continue gas injection. Injection would be discontinued and 
gas sales would begin to deplete the reservoir to about 
1300 psi or until compression was necessary to charge the sales 
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gas line. Production past this point was not considered part 
of this study. 

It was assumed that the wet gas produced would be in 
equilibrium at the pressure and temperature produced in the 
same proportion as that indicated by the equilibrium cell 
tests. In the equilibrium cell the recombined gas sample was 
allowed to expand at constant composition with each pressure 
decrement and then gas was extracted at constant pressure, 
thus simulating reservoir withdrawals. The material balance 
calculations resulting from these tests are believed to depict 
more realistically the actual reservoir performance than calcu- 
lation of the material balance by means of equilibrium constants. 

Based on a permeability variation of .30 as defined by 
Standing^ the dry gas cut at the producing wells was estimated 
from the first breakthrough at a producing well until the eco- 
nomic limit of injection was reached. 

Extraneous gas is available in excess of present avail- 
able market in a reservoir below the zone studied and conse- 
quantly at a higher pressure. This gas would be used as make 
up gas in the pressure maintenance case to compensate for 
losses due to fuel consumption and condensate production. 

After the economic limit is reached it is assumed that 
gas sales will begin again at the present rate and that pre- 
vailing prices would still be applicable. 

Since the field is just beginning to enter the retro- 
grade condensation phase considerably more condensate will be 
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produced by pressure maintenance if this procedure can be 
started at once. The greater the decline in pressure before 
injection commences the greater will be the losses in the 
reservoir and the higher will be the value of extraneous gas 
inj ec t ed. 

Because of the loss in gas sales during the period of 
pressure maintenance this procedure does not appear economi- 
cally attractive unless market availability for present gas 
sales can be filled from another source of gas. The ultimate 
recovery of condensate through pressure maintenance is excel- 
lent from a conservation point of view since over 9 0% of the 
available condensable hydrocarbon is expected to be recovered 
through the combination of dry gas injection to maintain pres- 
sure; to be followed by reservoir blowdown to depletion. 

C . Cycling with Pressure Depletion 

Return to the reservoir of available gas from the high 
pressure separator constitutes normal cycling. Since a net 
voidage situation exists as a result of shrinkage and fuel 
losses from the produced gas, normal cycling results in only 
partial pressure maintenance. However, the pressure reduc- 
tion would be little more than 407 o during the full economic 
cycling life. 

As in the case studied for full pressure maintenance, 
two wells would be used for injection of residue gas and four 
wells would be used for production of the wet gas effluent as 
indicated in the sweep model study. (See Figure 8.) 
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In this analysis gas sales are continued at the present 
rate which is dictated by the available gas market. Provision 
is made for compressors adequate to compress 40 million stand- 
ard cubic feet of residue gas per day for injection. All wet 
gas produced is processed through the low temperature units on 
offshore platforms. The sales gas and all the condensate pro- 
duced is charged to the sales gas line provided by the gas 
company , 

The amount of condensate produced will be governed by 
the equilibrium saturation at the prevailing producing bottom 
hole pressure until breakthrough of detectable dry gas. From 
this point onward the increasing dry gas cut will be governed 
by the permeability stratification as well as by the areal 
sweep pattern. As the cycling reaches its economic limit, the 
compressors will be shut down and sales gas production will 
continue to deplete the field. Production below 1300 psi was 
not considered in this study, since compression would be re- 
quired to charge the sales gas line which is maintained in 
excess of 1000 psi. 

It is predicted that normal cycling to .67 reservoir- 
gas cut followed by pressure blowdown will result in the recov- 
ery of 26,876,000 stock tank barrels of condensate or 67% of 
the available condensable hydrocarbons. Although the total 
amount of condensate produced through normal cycling is con- 
siderably less than for the pressure maintenance case, the 
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present worth of the income from early gas sales tended to 
make normal cycling more attractive economically. The ad- 
vantage of normal gas cycling, then, is continued income from 
the available market for gas while producing more condensate 
than would have been produced through normal pressure deple- 
tion and by producing this condensate earlier in the life of 
the field and in excess of available market for residue gas. 



CHAPTER V 



CONCLUSIONS 

A . Comparison of Ultimate Recoveries 

In a condensate reservoir, income is derived from the 
sales, not only of the residual gas produced, but also the 
liquids produced from the reservoir gas by mechanical sepa- 
ration or other type of plant processing. The ultimate 
quantity of these fluids produced is of importance as well 
as the rate of this production. 

The original gas in place is the same in any case. 

The ultimate produced quantity of this gas, however, is 
governed by the amount of fuel required in its production, 
or in the production of its condensate, and also by the reser- 
voir pressure at abandonment. The time at which this gas will 
be produced is dependent, of course, on the rate at which it 
is produced which in turn is dependent on the available mar- 
keting possibilities. In each of the alternatives studied the 
ultimate quantity of gas recovered was essentially the same, 
while the timing of the recovery of the gas was dependent on 
the method of field exploitation. 

The amount of condensate recovered is dependent on the 
pressure within the reservoir while being produced, since the 
retrograde properties of the reservoir fluid causes much of 
the heavier liquid hydrocarbon to remain within the reservoir 



45 



46 



rock. By producing the reservoir at a pressure at which retro- 
grade condensation will not occur or where it approaches some 
minimum value the quantity of condensate ultimately recovered 
is increased. 

The rate at which condensate is produced is dependent 
upon the rate of gas production unless some extraordinary 
means is provided to extract the condensate in advance of gas 
sales. This is the case in most gas reservoirs since pressure 
depletion without gas cycling is generally the method used for 
recovery. 

By cycling a part or all of the gas produced back into 
the reservoir not only can the condensate be extracted in 
advance of a gas market, but also the pressure can be main- 
tained at a level where little or no retrograde loss will 
occur in the reservoir. In the pressure maintenance case 
the rate of condensate recovery is limited only by the size 
of compressor plant installed. In the normal cycling case 
the rate is determined by the sum of the available gas market 
coupled with the capacity of the compressors. Under pressure 
maintenance essentially all, 95%, of the condensate is ulti- 
mately recovered and through normal cycling 67% is recovered. 
When compared with recovery of condensate by pressure deple- 
tion of only 3 8 % of the liquid, it is apparent that cycling 
is an important conservation measure. 

A comparison of ultimate recovery for both gas and 
condensate is presented in Table I. A plot of predicted 



cumulative condensate recovery versus time is presented as 
Figure 9 . 

B . Economic Analysis of Alternatives 

Since the quantity of gas ultimately to be recovered is 
essentially the same by any method of exploi tation ; the only 
factor in gas production contributing to its present value, 
and hence important in an economic analysis, is the time at 
which this gas is produced. In order to obtain maximum value 
from this valuable resource, the gas should be produced in as 
great a quantity as can be marketed. In the case of pressure 
depletion and in the case for normal gas cycling, gas was sold 
at the maximum rate limited only by the market obtainable. In 
the case for pressure maintenance the sale of gas was postponed 
for thirty years and this fact contributed to the reduced pres- 
ent worth for this alternative. 

Through pressure depletion 15,430,000 barrels of con- 
densate would be produced which is 38.6 per cent of the liquid 
originally in the reservoir. Should cycling while continuing 
to market gas be adopted, 26,876,000 million barrels of con- 
densate would be recovered. This is possible because increased 
condensate recovery while the reservoir pressure is still rela- 
tively high is made possible by the larger quantity of gas 
passing through the separators. Under the pressure maintenance 
program nearly all of the condensate, 38,090,000 barrels, can 
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Figure 9 
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be expected to be recovered, since the retrograde condensation 
caused by drop in reservoir pressure will not occur. 

A summary of the economic comparison of the three alter- 
native studies is included as Table I. From this summary it 
can be seen that the normal cycling program has the greatest 
present worth to an investor. This results from maximizing 
the value of gas in place by making use of the available mar- 
ket and combining this income with increased condensate recov- 
ery through retrieval at elevated pressure. Earlier retrieval 
of this condensate further contributes to its present net worth. 

Pressure maintnenace results in the maximum ultimate 
recovery of both the residue gas and the condensate. The post- 
ponement of gas sales until the economic limit of the cycling 
of the injection gas is reached, however, makes this procedure 
unattractive economically. Only if there was no market for 
gas at the present time, or if the available market could be 
satisfied by gas from another pool, could investment in pres- 
sure maintenance prove feasible. In this event, pressure 
maintenance would be the logical choice since the gas would 
have no present value for any of the alternatives and extra- 
neous gas injection to maintain pressure would result in 
valuable recovery of condensate not otherwise obtainable. 

The pressure depletion of the reservoir proves to be the 
least attractive program both from the conservationist viewpoint 
of quantity of hydrocarbons recovered and from the investor's 
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TABLE I 

ECONOMIC SUMMARY 





Pressure 

Depletion 


Pressure 
Ma i n t e nanc e 


Normal 

Cycling 


Pro j ec t Life 


50 year's 


80 years 


50 years 


Condensate Production, 
M- Bar r e 1 s 


15,430 


38,090 


26,876 


Total Sales Gas Production, 
MMMSCF 


492.7 


492.7 


492. 7 


Ultimate Condensate 

Recovery, % of original 


38. 6 


95.0 


67.2 


Gross Working Interest 
Income, M- $ 


114,934 


175,367 


142 , 860 


Operating Expense, M- $ 


4, 100 


16,367 


11,760 


Operating Income, M- $ 


110,834 


159 , 000 


131, 100 


Capital Investment, M- $ 


— 


4, 000 


5,000 


Total Cash Flow, M- $ 


110,834 


155,000 


126, 100 


Present Value, M- $ 


38,750 


30,000 


47,784 



5 1 



approach to maximizing profit from the invested dollar. From 
market data available pertaining to recovery from this field, 
some steps toward increasing condensate recovery by gas cycling 
should be taken at the earliest possible time. Since the res- 
ervoir is entering the retrograde phase now, time is of the 
essence if losses from this phenomenon are to be prevented. 

Graphs of the Annual and Cumulative Cash Plans versus 
time are shown for the three alternatives studied in Figures 10, 
11, and 12 . 

Capital expenditures required for gas cycling are given 



in Table VIII, page 72. 
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Figure 10 
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Figure 11 



Annual and Cumulative Cash Flow 
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Figure 12 



Annual and Cumulative Cash Flow 

Normal Gas Cycling with 
Pressure Depletion 
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CHAPTER VI 



RECOMMENDATIONS 

Based on the conclusions determined by the investigation 
of the supplementary condensate recovery programs studied, it 
is recommended that the following course of action be under- 
taken by the field operators: 

1. Organize an active operators committee for the 
express purpose of working toward the unitization of the 
reservoir. 

2. Charge the technical subcommittee with the respon- 
sibility of a detailed engineering design for a successful 
normal cycling program. 

3. Charge the appropriate subcommittee with the prepa- 
ration of unitization and operating agreements. 

4. Dedicate the manpower necessary to complete the 
engineering design at the earliest possible time to gain the 
advantage of the highest reservoir pressure attainable during 
recovery. 

5. Commence normal cycling of the residue gas to in- 
crease condensate recovery at the earliest practicable date. 
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DEFINITION OF TERMS 



We Wa ter influx 

C e Effective water compressibility (taking into account 

rock expansion into pore space because of geostatic 
compression of the reservoir rock 

Equivalent field radius 

Gp W Cumulative wet gas produced (SCF) 

G Initial gas in place (SCF) 

K Permeability in darcies 

Porosity 

Viscosity of water 
h e Effective pay thickness 

e v van Everdingen and Hurst water influx constant 

Qsc Flow rate for well tests 

Tf Temperature of formation (degrees Rankine) 

(j. g Viscosity of gas 

m Slope of drawdown curve (psi/cycle) 

P Reservoir pressure (psia) 

Bg Gas formation volume factor 

Z Gas deviation or compressibility factor 

i Subscript denoting initial value 

AP Reservoir pressure decrement 

Q r Infinite aquifer values of dimensionless water influx 

for values of dimensionless time, t ^ 

tj) Dimensionless time 

A Water influx constant 

T Real time 

Vg Bulk volume of the reservoir 
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TABLE II 

CALCULATION OF INITIAL GAS IN PLACE BY 
VOLUMETRIC METHOD 



Area 


Planime ter 
(sq . in. ) 


Area 
(acres ) 


Rat io 
o f 

Areas 


Interval 

ft. 


Equation 


AV 

acre-ft. 


A o 


21.60 


1983 . 46 










A 1 


15. 15 


1391. 18 


. 7 0 


20 


Trap 


33746 


A 2 


12.65 


1161.61 


. 83 


20 


Trap 


25528 


A 3 


10.23 


939.39 


. 81 


20 


Trap 


21010 


*4 


8.44 


775.02 


. 83 


20 


Trap 


17 144 


A 5 


6 .80 


624.42 


.80 


20 


Trap 


13994 


\ 


5 .35 


491.27 


. 79 


20 


Trap 


11157 


*7 


. 52 


47.75 


. 096 


20 


Pyr . 


4613 


A 8 


. 14 


12 . 86 


. 27 


20 


Pyr . 


569 


*9 


0 


0 


. 00 


15 


Pyr . 


64 



1 in^ = 91.827 acres 



h 

Pyramidal formula: AV R ~ ( A +A .-+-/A A ~ 

fi 3nn+lnn+l 

Trapezoidal formula: AV- = -x(A + A t1 ) 

B 2 n n+1 

V R = 127,825 acre-ft 
G = 7758 x V B x 4 X (1 - S w ) x B g 

G = 7758(.28)(.8)(1862.2)(127,825) = 415 x 10 9 SCF 
<t> = . 28 

S w = .20 

B g = .000537 bbl/cu. ft. or 1862.2 SCF/bbl 

Calculation of Flow Capacity 

(0.819)(Q gc )(T)(Z)(n g ) 

Kh e (tn)(P) 



60 



TABLE 1 1 - - Co n t i nu e d 






e 



(0.819)(9400)(700)(1.327)(. 034) 
(34.5 ) (8000) 



. 880 darcy- f t . 



K 



. 880 
17 . 5 



. 050 dare ies 



TABLE III 
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MONTHLY PAST PRODUCTION 





MCF GAS 


BBLS DISTILLATE 


BBLS WATER 


1/61 


76,227 


3,733 


8 


2/61 


250, 783 


15,098 


368 


3/61 


211,234 


16,196 


399 


4/61 


409,576 


27,841 


130 


5/61 


538,319 


42,100 


86 


6/61 


543,289 


42,403 


282 


7/61 


388,364 


28,872 


143 


8/61 


263, 335 


18,530 


35 


9/61 


217,357 


15,133 


84 


10/61 


433,201 


30,583 


138 


11/61 


553,913 


39,992 


300 


12/61 


729,040 


54,549 


217 


1961 


4,614,638 


335,030 


2,190 


1/62 


750,948 


54,465 


279 


2/62 


696,327 


52,138 


112 


3/62 


758,887 


57,031 


465 


4/62 


741,272 


56,523 


420 


5/62 


756,234 


55,536 


434 


6/62 


618,763 


45,731 


322 


7/62 


675,888 


52,441 


319 


8/62 


696,371 


47,460 


465 


9/62 


721,975 


49,439 


348 


10/62 


712, 711 


49,682 


496 


11/62 


787,513 


60,470 


420 


12/62 


788,886 


60,725 


527 


1962 


8,705,775 


641,641 


4,607 


1/63 


723,591 


61,729 


465 


2/63 


626,203 


49,738 


28 


3/63 


697,915 


51,312 


31 


4/63 


610,953 


49,704 


300 


5/63 


514,623 


40,434 


243 


6/63 


358,040 


31,177 


108 


7/63 


625,669 


47,156 


280 


8/63 


620,714 


49,127 


424 


9/63 


665,982 


45,531 


— 


10/63 


666,848 


46,491 


— 


11/63 


686,364 


43,249 


— 


12/63 


691,439 


49,899 


— 


1963 


7,470,341 


565,547 


1,879 


TOTAL 


20,790,754 


1,542,218 


8,676 
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TABLE IV 



PRESSURE DRAWDOWN-BUILDUP TEST 
Of f shore , Louisiana 



Pressure Drawdown 



Date Time 


Elapsed 

Time 

Days 


Bottom Hole 
Press. Psig. 


Remarks 


3-24-64 9:00 A.M. 






Wound 72 hour clock 


9:10 






Engaged clock 6c stylus 


9:15 






Pressured lubricator 


9:35 




8236 


Bomb placed at 11,745' 


10:45 


0 


8236 


Started opening well 


10:52-1/2 


.0052 


7943 


Well flowing 9400 MCF 


11:00 


.0104 


7892 




11:07-1/2 


.0156 


7884 




11:15 


.0208 


7884 




11:30 


.0313 


7876 




11:45 


.0417 


7876 




12:00 noon 


.0521 


7867 




12:30 P.M. 


.0729 


7867 




1:00 


.0938 


7859 




1:30 


.1146 


7867 




2:00 


.1354 


7867 




2:30 


.1562 


7867 




3:00 


.1771 


7867 




3:30 


.1979 


7867 




4:00 


.2187 


7867 




4:30 


.2396 


7867 




5:00 


.2604 


7867 




5:30 


. 2812 


7867 




6:00 


.3021 


7867 




6:30 


. 3229 


7859 




7:00 


.3437 


7859 




7:30 


.3646 


7859 




8:00 


.3854 


7859 




8:30 


.4062 


7859 




9:00 


.4271 


7859 




9:30 


.4479 


7859 




10:00 


.4687 


7851 




10:30 


.4896 


7851 




11 :00 


.5104 


7851 




11:30 


.5312 


7851 




12:00 midnt. 


.5521 


7842 
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TABLE IV- -Continued 



Pressure Drawdown 



* 

Date Time 



Elapsed 

Time 

Days 



Bottom Hole 
Press. Psig. 



Remarks 



3-25-64 12:30 A.M. 

1:00 
1:30 
2:00 
2:30 
3:00 
3:30 
4:00 
4:30 
5:00 
5:30 
6:00 
6:30 
7:00 
7:30 
8:00 
8:30 
9:00 
9:30 
10:00 
10:30 
11:00 
11:30 

12:00 noon 
12:30 P.M. 
1:00 
1:30 
2:00 
2:30 
3:00 
3:30 
4:00 
4:30 
5:00 
5:30 
6:00 
6:30 
7:00 
7:30 
8:00 
8:30 



.5729 


7842 


.5937 


7843 


.6146 


7825 


.6354 


7825 


.6562 


7817 


.6771 


7813 


.6979 


7809 


.7187 


7809 


. 7396 


7800 


.7604 


780Q. 


.7812 


7800 


.8021 


7792 


.8229 


7792 


.8437 


7792 


.8646 


7788 


.8854 


7788 


.9062 


7788 


.9271 


7784 


.9479 


7784 


.9687 


7784 


.9896 


7784 


1.0100 


7780 


1.0310 


7780 


1.0520 


7780 


1.0730 


7780 


1.0940 


7776 


1.1150 


7776 


1.1350 


7776 


1.1560 


7771 


1.1770 


7771 


1.1980 


7767 


1.2190 


7767 


1 . 2400 


7767 


1. 2600 


7763 


1.2810 


7759 


1.3020 


7759 


1.3230 


7755 


1 . 3440 


7755 


1.3650 


7751 


1.3850 


7751 


1.4060 


7746 
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TABLE IV- -Continued 



Pressure Drawdown 



Date 


Time 


Elapsed 

Time 

Days 


Bottom Hole 
Press. Psig. 


Remarks 


3-25-64 


9:00 P.M. 


1.4270 


7746 




(Cont . ) 


9:30 


1.4480 


7742 






10:00 


1.4690 


7742 






10:30 


1.4900 


7738 






11:00 


1.5100 


7738 






11:30 


1.5310 


7734 






12:00 midnt. 


1.5520 


7734 




3-26-64 


12:30 A.M. 


1.5730 


7734 






1:00 


1.5940 


7734 






1:30 


1.6150 


7734 






2:00 


1.6350 


7734 






2:30 


1.6560 


7730 






3:00 


1.6770 


7730 






3:30 


1.6980 


7725 






4:00 


1.7190 


7725 






4:30 


1.7400 


7725 






5:00 


1.7600 


7721 






5:30 


1.7810 


7721 






6:00 


1.8020 


7717 






6:30 


1.8230 


7717 






7:30 


1.8650 


7717 






8:30 


1.9060 


7709 






9:30 


1.9480 


7709 






10:30 


1.9900 


7700 






11:30 


2.0310 


7692 






12:30 P.M. 


2.0520 


7688 






1:30 


2.0730 


7683 






2:30 


2. 1560 


7680 






3:00 


2.1770 


7680 






4:00 


2.2190 


7675 






5:00 


2.2600 


7675 






6:00 


2.3020 


7671 






7:00 


2.3440 


7671 






8:00 


2. 3850 


7671 






9:00 


2.4270 


7671 


60 hrs. Shut well in. 
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TABLE IV- -Continued 



Pressure Buildup 



Time 




Elapsed 

Hours 


Time 

Days 


Surface DWT 
Press. Psig. 


Bottom Hoi 
Press. Psi, 


9:00 


P.M. 


0 


0 


4745 


7671 


9:01 




- 


.0007 


6480 


-- 


9:02 




- 


.0014 


6510 


-- 


9:03 




- 


.0021 


6530 


-- 


9:04 




- 


.0028 


6535 




9:05 




- 


.0034 


6542 




9:07- 


■1/2 


.125 


.0052 


6547 


7960 


9:10 




.167 


.0069 


6551 




9:15 




.250 


.0104 


6552 


7977 


9:22- 


-1/2 


. 370 


.0156 


6556 


6985 


9:30 




.500 


.0208 


6555 


7989 


9:37- 


-1/2 


.620 


.0260 


6555 


7994 


9:45 




.750 


.0312 


6554 


7994 


10:00 




1.000 


.0417 


6554 


7994 


10:15 




1.250 


.0521 


6552 


7994 


10:30 




1.500 


.0625 


6551 


7994 


10:45 




1.750 


.0729 


6550 


7994 


11:00 




2.000 


.0833 


6548 


7994 


11:30 




2.500 


. 1042 


6548 


7994 


12:00 


A.M. 


3.000 


.1250 


6548 


7994 


12:30 




3.500 


.1458 


6548 


7994 


1:00 




4.000 


. 1667 


6549 


7994 


1:30 




4.500 


. 1875 


6550 


7994 


2:00 




5.000 


.2083 


6552 


7994 


2:30 




5.500 


. 2292 


6555 


7989 


3:00 




6.000 


. 2400 


6562 


7985 


3:30 




6.500 


.2708 


6566 


7977 


4:00 




7.000 


.2917 


6570 


7977 


4:30 




7.500 


.3125 


6575 


7977 


5:00 




8.000 


.3333 


6578 


7977 


5:30 




8.500 


.3542 


6579 


7973 


6:00 




9.000 


.3750 


6580 


7968 


6:30 




9.500 


.3958 


6583 


7968 


7:00 




10.000 


.4167 


6585 


7968 


7:30 




10.500 


.4375 


6586 


7968 


8:00 




11.000 


.4583 


6586 


7968 


8:30 




11.500 


.4792 


6586 


7968 


9:00 




12.000 0 


.5000 


6586 


7968 



Remarks 

Shut well in. 



Started out 
of hole. 
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TABLE IV- -Continued 



Static Pressures 
15,000# Element 



Depth 

0 

10.500 

11,000 

11.500 
11,780 
11,835 



Bottom Hole 
Press. Psig. 

6627 

7993 

8018 

8068 



AP 



25 

50 



Gradient (Psi/ft.) 



.050 

.100 



8102 



34 



. 121 



Year 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 



TABLE V 



PRESENT WORTH 

STRAIGHT DEPLETION OF PRESSURE 
AT CURRENT RATE 



Condensate Wet Gas Gross 
Bbls MMSCF Income 



617490 


9346 


3052287 


585247 


9346 


2973627 


547395 


9346 


2881280 


511787 


9346 


2794409 


482908 


9346 


2623995 


456178 


9346 


2658742 


429448 


9346 


2593532 


409448 


9346 


2543778 


389447 


9346 


2495943 


272718 


9346 


2455130 


356082 


9346 


2414540 


342718 


9346 


2381940 


333745 


9346 


2360050 


318231 


9346 


2322200 


311502 


9346 


2305784 


300380 


9346 


2278650 


295987 


9346 


2267933 


289259 


9346 


2241420 


282623 


9346 


2235330 


278137 


9346 


2224385 


273744 


9346 


2213668 


271501 


9346 


2208195 


267015 


9346 


2197250 


262623 


9346 


2186533 


260380 


9346 


2181064 


258136 


9346 


2175589 


255893 


9346 


2170117 


255893 


9346 


2170117 


253650 


9346 


2164646 


251407 


9346 


2159172 



Oper . 
Costs 


Undeferred 
Net Income 


Deferred 
Net Profit 


124439 


2927848 


2842940 


120005 


2853622 


2616770 


114801 


2766479 


2390240 


109904 


2684505 


2190556 


105934 


2618061 


2013290 


102258 


2556484 


1856000 


98593 


2494949 


1790040 


95834 


2447944 


1581370 


93084 


2402859 


1465640 


90784 


2364346 


1359500 


88494 


2326046 


1263040 


86658 


2295282 


1175185 


85424 


2274626 


1098645 


83291 


2238909 


1020940 


82364 


2223420 


953850 


80836 


2197814 


890150 


80234 


2187699 


835700 


79307 


2172213 


784170 


78394 


2156936 


733360 


77778 


2146607 


689060 


77174 


2136494 


647360 


76864 


2131331 


609460 


76249 


2121001 


572670 


75644 


2110889 


563165 


75334 


2104730 


505375 


75028 


2100561 


474725 


74719 


2094398 


448415 


74719 


2095398 


423875 


74411 


2090235 


397145 


74102 


2085070 


375310 



31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 
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TABLE V--Continued 



Condensate 

Bbls 


Wet Gas 
MMSCF 


Gross 

Income 


Oper . 
Costs 


Undeferred 
Net Income 


Deferred 
Net Profit 


250380 


9346 


2146667 


73961 


2082706 


347810 


249258 


9346 


2153930 


73807 


2080123 


330740 


248136 


9346 


2151190 


13653 


2077537 


313710 


244770 


9346 


2142980 


73190 


2069790 


293910 


244770 


9346 


2142980 


73190 


2069790 


277351 


244770 


9346 


2142980 


73190 


2069790 


262860 


244770 


9346 


2142980 


73190 


2069790 


246300 


244770 


9346 


2142980 


73190 


2069790 


233890 


244770 


9346 


2142980 


73190 


2069790 


219400 


244770 


9346 


2142980 


73190 


2069790 


206979 


244770 


9346 


2142980 


73190 


2069790 


195180 


244770 


9346 


2142980 


73190 


2069790 


184830 


244770 


9346 


2142980 


73190 


2069790 


173860 


244770 


9346 


2142980 


73190 


2069790 


164340 


244770 


9346 


2142980 


73190 


2069790 


154820 


244770 


9346 


2142980 


73190 


2069790 


146334 


244770 


9346 


2142980 


73190 


2069790 


140124 


244770 


9346 


2142980 


73190 


2069790 


130600 


244770 


9346 


2142980 


73190 


2069790 


123150 


244770 


9346 


2142982 


73190 


2069790 


115910 


15429856 


467300 








38749249 
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TABLE VI 
PRESENT WORTH 

PRESSURE MAINTENANCE WITH GAS CYCLING 
40 MMSCF PER DAY 



Years 


Condensate 

Bbls 


Wet Gas 
MMSCF 


Gross 

Income 


Oper . 
Costs 


Undeferred 
Net Income 


Deferred 
Net Profit 


i 


1000000 




2439650 


441500 


1998150 


(-2059800) 


2 


1000000 




2439650 


441500 


1998150 


(-227497) 


3 


1000000 




2439650 


441500 


1998150 


1498900 


4 


1000000 




2439650 


441500 


1998150 


1630490 


5 


1000000 




2439650 


441500 


1998150 


1536600 


6 


1000000 




2439650 


441500 


1998150 


1450660 


7 


1000000 




2439650 


441500 


1998150 


1368700 


8 


1000000 




2439650 


441500 


1998150 


1290800 


9 


1000000 




2439650 


441500 


1998150 


1218900 


10 


1000000 




2439650 


441500 


1998150 


1149000 


11 


1000000 




2439650 


441500 


1998150 


1085000 


12 


1000000 




2439650 


441500 


1998150 


1023000 


13 


1000000 




2439650 


441500 


1998150 


965100 


14 


1000000 




2439650 


441500 


1998150 


911160 


15 


1000000 




2439650 


441500 


1998150 


857200 


16 


1000000 




2439650 


441500 


1998150 


809250 


17 


1000000 




2439650 


441500 


1998150 


763300 


18 


1000000 




2439650 


441500 


1998150 


721300 


19 


1000000 




2439650 


441500 


1998150 


679370 


20 


920000 




2244500 


430500 


1814000 


641400 


21 


890000 




2171300 


426375 


1744925 


605440 


22 


860000 




2098100 


422250 


1675850 


571500 


23 


830000 




2024900 


418125 


1606775 


539500 


24 


810000 




1976100 


415375 


1560725 


507500 


25 


780000 




1902900 


411250 


1491650 


479560 


26 


750000 




1829700 


407125 


1422575 


451580 


27 


730000 




1780900 


404375 


1376525 


427600 


28 


710000 




1732200 


401625 


1330575 


403630 


29 


690000 




1683400 


398875 


1284525 


379650 


30 


670000 




1634600 


396125 


1238475 


359670 


IBTOTAL 


27640000 










24325760 


RESERVOIR BLOWDOWN 
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10450000 


482.3 


105625267 


3446150 


102179117 


5687293 


'OTAL 


38090000 


482.3 








30013053 



Year 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 
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TABLE VII 
PRESENT WORTH 

GAS CYCLING WITH PRESSURE DEPLETION 
GAS SALES AT CURRENT RATE 



Condensate 

Bbls 


Wet Gas 
MMSCF 


Gross 

Income 


1605475 


9345.6 


5462600 


1521642 


9345.6 


5258100 


1423230 


9345.6 


5018000 


1330646 


9345.6 


4792100 


1255560 


9345.6 


4608950 


1186062 


9345.6 


4439400 


1116565 


9345.6 


4269860 


1064565 


9345.6 


4143000 


1012562 


9345.6 


4016130 


969067 


9345.6 


3910000 


925813 


9345.6 


3804500 


891067 


9345.6 


3719700 


867737 


9345.6 


3662800 


744660 


9345.6 


3362540 


688420 


9345.6 


3225330 


632300 


9345.6 


3088420 


588720 


9345.6 


2982100 


545253 


9345.6 


2876060 


510700 


9345.6 


2791760 


477280 


9345.6 


2710244 


273744 


9345.6 


2213670 


271500 


9345.6 


2208200 


267015 


9345.6 


2197250 


262623 


9345.6 


2186540 


260380 


9345.6 


2181060 


258136 


9345.6 


2175600 


255893 


9345.6 


2170120 


255893 


9345.6 


2170120 


253650 


9345.6 


2164650 


241407 


9345.6 


2149170 



Oper . 
Cos ts 


Undeferred 
Net Income 


Deferred 
Net Profit 


564285 


4898315 


(-243735) 


552760 


4705340 


4071065 


539230 


4478770 


3869660 


526500 


4265600 


3480730 


516175 


4092775 


3147350 


506620 


3932780 


2855200 


497065 


3772795 


2584360 


389913 


3653086 


2359900 


482765 


3533365 


2155350 


476785 


3433215 


1974100 


470835 


3333665 


1810180 


466055 


3253645 


1665870 


462845 


3199955 


1545600 


445925 


2916615 


1330000 


439193 


2787137 


1195700 


430475 


2657945 


1076467 


424485 


2557615 


977000 


418535 


2456525 


887170 


413755 


2378005 


808521 


409165 


2301079 


738810 


77175 


2136495 


647360 


76865 


1231335 


609560 


76265 


2120985 


572670 


75645 


2110895 


536170 


75335 


2105725 


505375 


75025 


2100575 


474730 


74720 


2095400 


448415 


74720 


2095400 


523270 


74410 


2090240 


397145 


74105 


2085065 


375310 



TABLE VII--Continued 



Year 


Condensa te 
Bbls 


Wet Gas 
MMSCF 


Gross 

Income 


Oper . 
Costs 


Undeferred 
Net Income 


Deferred 
Net Profit 


31 


250380 


9345.6 


2156670 


73965 


2082705 


347810 


32 


249248 


9345.6 


2153930 


73805 


2080125 


330740 


33 


248136 


9345.6 


2151190 


73665 


2077535 


313700 


34 


244770 


9345.6 


2143000 


73190 


2069810 


293910 


35 


244770 


9345.6 


2143000 


73190 


2069810 


277350 


36 


244770 


9345.6 


2143000 


73190 


2069810 


262870 


37 


244770 


9345.6 


2143000 


73190 


2069810 


246300 


38 


244770 


9345.6 


2143000 


73190 


2069810 


233900 


39 


244770 


9345.6 


2143000 


73190 


2069810 


219400 


40 


244770 


9345.6 


2143000 


73190 


2069810 


207000 


41 


244770 


9345.6 


2143000 


73190 


2069810 


195180 


42 


244770 


9345.6 


2143000 


73190 


2069810 


184800 


43 


244770 


9345.6 


2143000 


73190 


2069810 


173860 


44 


244770 


9345.6 


2143000 


73190 


2069810 


164340 


45 


244770 


9345.6 


2143000 


73190 


2069810 


154820 


46 


244770 


9345.6 


2143000 


73190 


2069810 


146340 


47 


244770 


9345.6 


2143000 


73190 


2069810 


140130 


48 


244770 


9345.6 


2143000 


73190 


2069810 


130600 


49 


244770 


9345.6 


2143000 


73190 


2069810 


123150 


50 


244770 


9345.6 


2143000 


73190 


2069810 


115900 



TOTAL 26876429 467280 47784138 
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TABLE VIII 

CAPITAL EXPENDITURES 



Pressure Maintenance - Cycling without gas sales 



Re comp 1 e t e 9 
New Well 11 
Recomplete 10 
N ew Well 5 
Recomp 1 e t e 4 
Well 6 



Platform 

Fabrication 

Installation 



Compressors 

Fabrication 
Ins tallation 



Pipe 

High pressure line 
Sales gas line 



Separation Facilities 



$ 70,000 

600, 000 
70, 000 
600, 000 
70, 000 
0 



$1,410,000 



320, 000 
8 0 , 000 

$ 400,000 



1,116, 000 

90, 000 

$1,206, 000 



25, 000 

100, 000 

$ 125,000 
$ 971,500 
$4, 000, 000 



Additional Plant Required for Normal Cycling 

New well $ 700, 000 

Sales gas line 300, 000 

$ 1 , 000,000 

Total for normal gas cycling $5 , 000, 00 0 



APPENDIX I 



MATERIAL BALANCE CALCULATION 



GAS RESERVE CALCULATION PROGRAM 
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800 

101 



200 

1000 



201 

1010 

c 



202 

1020 

C 



203 

C 



204 

1030 

102 



300 



PROGRAM GA5RES 

DIMENSION A ( 6 » 6 ) *B(6) *T (381) » Q ( 38 1 ) ,DELP(20) , P ( 20 ) * GP ( 2 0 ) , DP ( 20 
DIMENSION WP ( 20 ) ,GPW(20> »WI (20) »DELBG ( 20 ) 

COMMON A,B,TEMP,PC,TC,T,Q,K,DELP 
PRINT 800 
FORMAT ( 1H1 ) 

PRINT 101 

FORMAT ( 44H RESERVOIR COMPUTER PROJECT ENERGY BALANCE,//) 

DATA INPUT FOR SUBROUTINE ZCALC CONSTANTS 
DO 1000 J= 1 » 6 » 1 
DO 1000 1=1 , 6.1 
READ 200, A( I , J) 

FORMAT ( 13X.F14.il) 

CONTINUE 
DO 1010 1=1, 6,1 
READ 201 ,B( I ) 

FORMAT ( F14.ll ) 

CONTINUE 

DATA INPUT FOR SUBROUTINE QFUNCT TABLE 
DO 1020 1=1,381,3 

READ 202,T ( I ) ,Q( I ) ,T ( 1 + 1 ) ,Q( 1+1 ) ,T ( 1+2 ) ,Q( 1+2 ) 

FORMAT ( FI 2. 2, FI 2. 3, FI 2. 2, FI 2. 3, FI 2. 2, FI 2. 3 ) 

CONTINUE 

DATA INPUT FOR RESERVOIR PROPERTIES 

READ 203»AIWEL 

READ 203 »CONV 

READ 203, RI 

READ 203, DEPTH 

READ 203, TEMP 

READ 203, SGI 

READ 203,0 

READ 203, AH 

READ 203, SG 

READ 203, U 

READ 203, AK 

READ 203, TIME 

READ 203, SGR 

READ 203, SWC 

READ 203, PI 

READ 203, PC 

READ 203, TC 

READ 203, BW 

READ 203, CPI 

FORMAT ( FI 6 . 7 ) 

I NC=CP I 

DATA INPUT FOR RESERVOIR PRODUCTION 

DO 1030 1 = 1 , INC, 1 

READ 204, GP( I) , DP ( I ) , WP ( I ) 

FORMAT (3F20.0) 

CONTINUE 
PRINT 102 

FORMAT ( 1 9H RESERVOIR DATA FOR,/) 

DATA PRINT OUT FOR RESERVOIR PROPERTIES 
PRINT 300, INC 

FORMAT ( 45H INC, CURRENT NUMBER OF PRODUCTION I NCREMENTS= , 14 ) 
PRINT 301 , A I WEL 
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301 FORMAT ( 48H MAX. NUMBER OF PRODUCING WELLS AT ANY ONE TIME=»F6.1) 
PRINT 302, BW 

302 FORMAT ( 50H BW, FORMATION VOLUME FACTOR FOR WATER IN BBLS /STB = , F6 . 3 ) 
PRINT 303 » CONV 

303 FORMAT ( 44H CONV,SCF OF VAPOR PER BARREL OF D I ST I LL AT E= , F7 . 1 ) 

PRINT 304, RI 

304 FORMAT ( 48H RI, RADIUS TO INITIAL EDGE OF RESERVOIR IN FEET=,F8.1) 
PRINT 305, DEPTH 

305 FORMAT ( 34H DEPTH, DEPTH OF RESERVOIR IN FEET=,F8.1) 

PRINT 306, TEMP 

306 FORMAT ( 50H TEMP , TEMPERATURE OF RESERVOIR IN DEGREES RANK I NE= , F7 . 1 ) 
PRINT 307 , PC 

307 FORMAT ( 36H PC ,PSEUDOCR I T ICAL PRESSURE IN PSIA=,F6.1) 

PRINT 308 , TC 

308 FORMAT ( 50H TC ,PSEUDOCR I T ICAL TEMPERATURE IN DEGREES RANK I NE= » F6 • 1 ) 
PRINT 309, SGI 

309 FORMAT ( 40H SGI, INITIAL GAS SATURATION AS FRACT I ON= , F5 . 3 ) 

PRINT 310,0 

310 FORMAT ( 24H 0, POROSITY AS FRACT I ON= , F5 . 3 ) 

PRINT 311, AH 

311 FORMAT ( 43H AH, AVERAGE THICKNESS OF RESERVOIR IN FEET=,F6.1) 

PRINT 312, SG 

312 FORMAT ( 44H SG, SPECIFIC GRAVITY OF GAS RELATIVE TO AIR=,F7.4) 

PRINT 313, U 

313 FORMAT ( 46H U, VISCOSITY OF RESERVOIR WATER IN CENT I POI SE= , F6 . 3 ) 
PRINT 314, AK 

314 FORMAT ( 47H AK, WATER PERMEABILITY OF RESERVOIR IN DA RC I ES= , F7 . 4 ) 
PRINT 315, TIME 

315 FORMAT ( 33H TIME, TIME INTERVAL USED IN DAYS=,F7.2) 

PRINT 316, SGR 

316 FORMAT ( 41H SGR, RESIDUAL GAS SATURATION AS FRACT I ON = * F5 . 3 ) 

PRINT 317, SWC 

317 FORMAT ( 42H SWC, CONNATE WATER SATURATION AS FRACT I ON= , F 5 . 3 ! 

PRINT 802 

802 FORMAT ( 1H ,/) 

! DATA INPUT FOR RESERVOIR PRESSURE 

1 DO 1040 1=1 , INC, 1 
READ 205 ,P ( I ) 

205 FORMAT ( F10. 3 ) 

1040 CONTINUE 

: MAIN BODY OF PROGRAM WITH CALCULATIONS 

PRINT 103 

103 FORMAT ( 16H TIME PR ESSURE , 9X , 1HZ , 1 OX , 1 OHBG , BBL /SC F , 8X , 6HDEL BG,6X 
1,1 4HRES • WITH., BBL) 

DELP ( 1 ) =0.0 
XW=0.0 
YW=0 . 0 
XY=0 • 0 
Y Y=0 • 0 
WW=0.0 

C=0.00514590*TEMP 

DELTD=1858064.0*AK*TIME/ (U*0*RI*RI ) 

S= P ( 1 ) 

CALL ZCALC ( S , Z ) 

BGI=C*Z/S 
Y = 0 . 0 
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L = 0 
X = 0.0 

PRINT 400,L»S,Z,BGI >Y,X 

400 FORMAT ( I5,2X,F10.3,2X,E14.8,2X,E14.8,2X,E14.8,2X,E14.8) 

DO 2 K=2*INC»1 

L = K-1 

DELP(K)=P(K-1 )-P(K) 

CALL INFLUX (DELTD»WD) 

WI ( K ) =WD 
W=WD 
S=P ( K ) 

CALL ZC ALC ( S * Z ) 

BG=C#Z/S 

GPW ( K ) =GP ( K ) *BG+WP ( K ) *BW+DP ( K ) *CONV*BG 
X=GPW ( K ) 

Y= BG-BG I 
DELBG ( K ) =Y 
XW=X*W+XW 
X Y=X*Y+XY 
YW=Y*W+YW 
YY=Y*Y+YY 

ww=w*w+ww 

PRINT 400»L»S»Z»BG»Y »X 

2 CONTINUE 

G= ( ( XY*WW )-( XW*YW )) / I ( WW*YY )-( YW*YW) ) 

AA=(XW-G*YW)/WW 
PRINT 40 1 » G » A A 

401 FORMAT ( 1H *//*14H L.M.S. G > SCF= , E 14 . 8 , / , 23H WATER INFLUX CONSTANT = 
2 » E 14. 8 » / / ) 

PRINT 104 

104 FORMAT ( 1 6H TIME PRESSURE » 6X » 5HDEL P * 3X > 14HWATER INFL , BBL » 2X * 14HP 
3ERIODIC G CHK ,2X , 12HPERCENT DEV.) 

WE = 0 • 0 
L = 0 

PRINT 402»L*PI»DELP(1 ) ,WE 

402 FORMAT ( I5,2X,F10.3,2X,F10.3,2X>E14.8) 

DO 3 K= 2 » I NC ♦ 1 

L = K- 1 

WE=AA#W I (K) 

GC=(GPW(K)~WE)/DELBG(K) 

PCENT=( GC“G)*100.0/G 

PRINT 403 » L » P ( K ) > DELP (K ) »WE»GC»PCENT 

403 FORMAT ( I5>2X,F10.3»2X,F10.3,2X,E14.8,2X,E14.8,2X,F10.4) 

3 CONTINUE 
PRINT 801 

801 FORMAT ( 1H ,/////) 

PRINT 800 
CONTINUE 
GO TO 1 
END 
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SUBROUTINE ZCALC(S,Z) 

DIMENSION A ( 6 > 6 ) >B ( 6 ) » T ( 38 1 ) »Q ( 38 1 ) »DELP<20) • P C 20 ) > GP ( 2 0 ) , DP ( 2 0 ) 
DIMENSION PX ( 6 ) ,PY( 6 ) 

COMMON A»B»TEMP»PC»TC»T >Q > K. » DELP 
Y=( (2.0*TEMP/TC)-4.0)/1.9 
PY ( 1 ) =B ( 1 ) 

PY ( 2 ) =B( 2 ) *Y 

PY( 3 ) =B ( 3 ) *( { 3.0*Y*Y ) -1 .0 ) 

PY ( 4 ) =B ( 4 ) * ( ( 5.0*Y*Y) -3#0)*Y 
P Y ( 5 ) =B ( 5 ) * ( ( ( 35.0*Y*Y 1-30.0 ) *Y*Y+3.0 ) 

PY(6)=B(6)*((< (63.0*Y*Y)-70.0)*Y*Y)+15.0)*Y 
X=( (2.0*S/PC)-15.0)/14.8 
PX ( 1 ) =B( 1 ) 

PX ( 2 ) =B ( 2 ) *X 

PX ( 3 ) =B ( 3 ) * ( ( 3 • 0*X*X 1-1.0) 

PX(4)=B(4)*( ( 5.0*X*X )-3.0)*X 
PXf5)=B(5)*(({ (35.0*X*X ) -30.0) *X*X 1+3.0) 

PX(6)=B(6)*((( (63.0*X*X 1-70.0 )*X*X )+15.0)*X 
APP=0 . 0 
DO 1 M= 1 » 6 » 1 
DO 1 N=l»6»l 

1 APP=( A(N»M)*PX(M)*PY( N) )+APP 
Z = APP 
RETURN 
END 



81 



SUBROUTINE OFUNCT ( DELT ,QT ) 

DIMENSION A ( 6 ♦ 6 ) > B ( 6 ) » T ( 38 1 ) » Q ( 3 8 1 ) ,DELP(20) ,P(20)*GP(20),DP(20) 

COMMON A»B,TEMP»PC,TC»T ,Q»K»DELP 

TT=DELT 

DO 3 1=1*381*1 
I F ( TT-T ( I ) ) 1,2,3 

1 DQ1 = (Q(I-1) -Q ( 1-2) )/(T( I-l)-T( 1-2) ) 

DQ2 = ( (Q(I)-Q(I-1))/(T(I)-T(I-1) )-DQl ) / ( T ( I ) -T ( I -2 ) ) 

QT=Q( I — 2 ) + ( TT-T ( 1-2) )*DQ1+(TT-T( 1-2) )*(TT-T( I - 1 ) )*DQ2 
GO TO 4 

2 QT = Q ( I ) 

GO TO 4 

3 CONTINUE 

4 RETURN 
END 

END GASRES 
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SUBROUTINE INFLUX! DFLTD»WD) 

DIMENSION A ( 6 » 6 ) » B ( 6 ) » T ( 38 1 ) *0(381) >DELP(20) ,P(20)»GP(20),DP(20) 
COMMON A,B,TEMP,PC,TC,T,Q,K,DELP 
WD = 0 . 0 

DO 1 J= 1 » K ♦ 1 
T I = J- 1 

DELT=TI*DELTD 
CALL QFUNCT ( DELT »QT ) 

N=K+1-J ' 

WI=DELP(N)*QT 
WD=W I +WD 
1 CONTINUE 
RETURN * 

END 
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RESERVOIR COMPUTER PROJECT — “ENERGY BALANCF 



RESERVOIR DATA FOR 

INOCURRENT NUMBER OF PRODUCTION I NCREMENTS= 13 

MAX. NUMBER OF PRODUCING WELLS AT ANY ONE T I ME= 2.0 

BW ^FORMAT I ON VOLUME FACTOR FOR WATER IN BBLS/STB= 1.030 

CON V , SC F OF VAPOR PER BARREL OF D I ST I LLATE= 768.0 

RI, RADIUS TO INITIAL EDGE OF RESERVOIR IN FEET= 2190.0 

DEPTH »DEPTH OF RESERVOIR IN FEET= 11450.0 

TEMP, TEMPERATURE OF RESERVOIR IN DEGREES RANK I NE= 690.0 

PC,PSEUDOCRI TICAL PRESSURE IN PSIA= 662.3 

TC,PSEUDOCRI TICAL TEMPERATURE IN DEGREES RANK I NE= 387.5 

SGI, INITIAL GAS SATURATION AS FRACT I ON= .800 

0, POROSITY AS FRACT I 0N= .280 

AH, AVERAGE THICKNESS OF RESERVOIR IN FEET= 64. A 
SG, SPECIFIC GRAVITY OF GAS RELATIVE TO AIR= .7975 
U, VISCOSITY OF RESERVOIR WATER IN CENTIPOISE= .300 
AK, WATER PERMEABILITY OF RESERVOIR IN DARCIES= .0500 
TIME, TIME INTERVAL USED IN DAYS= 90.00 
SGR, RESIDUAL GAS SATURATION AS FRACT I ON= .300 
SWC»CONNATE WATFR SATURATION AS FRACTI0N= .150 



ME 


PRESSURE 


Z 




BG»BBL/SCF 


RES. WITH. 


, BBL 


0 


8838.000 


1. 32981926E 


00 


5.34255563E-04 




0 


1 


8818.000 


1 . 32803329E 


00 


5. 34748162E-04 


3 < 03008359E 


05 


2 


8764.000 


1 . 323220 17E 


00 


5.36093049E-04 


1. 14994894E 


06 


3 


8717.000 


1.3 1 90409 5E 


00 


5.3728 1226E-04 


1 . 645 49602 E 


06 


4 


8657.000 


1.3 1371808E 


00 


5.38821841 E-04 


2 . 62 736394E 


06 


5 


8582.000 


1. 30708189E 


00 


5.40785103E-04 


3. 89883102 E 


06 


6 


8509.000 


1.30063895E 


00 


5.42736044E-04 


5 . 128442 13 E 


06 


7 


8440.000 


1.29456183E 


00 


5.44616486E-04 


6 . 35037990E 


06 


8 


8364.000 


1.28788050E 


00 


5.46728833E-04 


7 . 69974695 E 


06 


9 


8289.000 


1 • 2 8 1 2978 1 E 


00 


5.48855951E-04 


8.91284897E 


06 


10 


8219.000 


1.27516200E 


00 


5.50879758E-04 


9 . 8 1497565E 


06 


11 


8155.000 


1.26955778E 


00 


5.52762967E-04 


1 .09665134E 


07 


12 


8100.000 


1 . 26474529E 


00 


5.54406719E-04 


1 .21921223E 


07 


M.S. 


G » SCF = 6 • 066 


24695E 11 











WATER INFLUX CONST ANT = 4 . 0706 80 1 #E 00 



ME 


PRESSURE 


DEL P 


WATER INFL,BBL 


PERIODIC G CHK 


1 


8818.000 


20.000 




0 


6. 15121449E 


11 


2 


8764.000 


54.000 


1 • 03 1 1 549*E 


03 


6.26388336E 


11 


3 


8717.000 


47.000 


4. 5347547*E 


03 


5.45345143E 


11 


4 


8657.000 


60.000 


9.5550091*E 


03 


5.77476619E 


11 


5 


8582.000 


75.000 


1 .6725031*E 


04 


5.99668013E 


11 


6 


8509.000 


73.000 


2.6552691*E 


04 


6.07865842E 


11 


7 


8440.000 


69.000 


3. 8609296 *E 


04 


6 . 16642836E 


11 


'8 


8364.000 


76.000 


5 . 2608605*E 


04 


6.21517492E 


11 


9 


8289.000 


75.000 


6. 8905530 *E 


04 


6.15172299E 


11 


10 


8219.000 


70.000 


8 . 7297406*E 


04 


5.95654287E 


11 


11 


8155.000 


64.000 


1 .0748285*E 


05 


5.98354920E 


11 


12 


8100.000 


55.000 


1 .2918882*E 


05 


6.1 144437 6E 


11 



APPENDIX II 



FLUID DATA 
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HYDROCARBON ANALYSES OF SEPARATOR PRODUCTS 
AND CALCULATED WELL STREAM 



Separator Liquid 


Separator Gas 


Well Stream 




GPM@ 


GPM@ 



Component 


Mol 7, 


Liq. Vol . ° 


\ Mol 7, 


15.025 PSIA 


Mol % 


15.025 PSIA 


Carbon Dioxide 








0.70 




0.67 




Nitrogen 


— 


— 


0.32 




0.30 




Methane 


21.59 


5.59 


93.34 




89.90 




Ethane 


4.13 


1.44 


3.53 




3.56 




Propane 


3.87 


1.95 


1.24 


0.349 


1.37 


0.385 


Iso-Butane 


2.05 


1.23 


0.35 


0.117 


0.43 


0.144 


N-Butane 


2.29 


1.32 


0.25 


0.081 


0.35 


0.113 


Iso-Pentane 


1.47 


0.99 


0.11 


0.041 


0.17 


0.063 


N-Pentane 


1.68 


1.12 


0.05 


0.018 


0.13 


0.048 


Hexanes 


4.82 


3.63 


0.07 


0.029 


0.30 


0.126 


Heptanes Plus 


58.10 


82.73 


0.04 


0.019 


2.82 


2.233 


Totals 


100.00 


100.00 


100.00 


0.654 


100.00 


3.112 


Calculated Specific Gravity (Air - 


1.00) 


0.606 




0.7975 



Properties of Heptanes Plus : 



API Gravity 
Specific Gravity 
Molecular Weight 



39 . 9° API @ 60°F 
0.8256 @ 60/60°F 
202 



Basis of Recombination: 



Separator Liquid per MMSCF Separator Gas = 66.07 Bbls 



PRESSURE-VOLUME RELATIONS OF RESERVOIR FLUID AT 240° F 
(CONSTANT COMPOSITION EXPANSION) 
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Pressure 

(PSIA) 


Relative 

Volume 

(V/Vsat) 


Specific 

Volume 

(Cu.Ft./Lb.) 


Deviation 

Factor 

(Z) 


8680 


0.9793 


0.04918 


1.3111 


8600 Original Reservoir 
Pressure 








8305 


0.9978 


0.05011 


1.2782 


8260 Dew Point Press. 


1.0000 


0.05022 


1.2741 


8138 


1.0070 


0.05057 




7825 


1.0255 


0.05150 




2555 


1.0442 


0.05244 




6930 


1.0910 


0.05479 




6415 


1.1382 


0.05716 




5990 


1.1856 


0.05954 




5610 


1.2332 


0.06193 




5295 


1.2808 


0.06432 




4765 


1.3767 


0.06914 




4168 


1.5207 


0.07637 




3495 


1.7555 


0.08816 




2675 


2. 2401 


0.11250 




2190 


2.7208 


0.13663 




1855 


3.2067 


0.16104 




1620 


3.6864 


0.18513 





EFFECT OF GAS-LIQUID RATIO UPON DEW POINT PRESSURE AT 240°F 



Dew Point Pressure 
(PS I A) 



8,090 

8,260 

8,360 



Gas-Liquid Ratio 
(SCF Separator Gas per 
Bbl Separator Liquid) 

20,000 

15,135 

10,000 
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RETROGRADE CONDENSATION DURING GAS DEPLETION AT 240°F 



Pressure 

(PSIA) 



Reservoir Liquid 
(Bbls/MMSCF of 
Reservoir Fluid) 



Retrograde 
Liquid Volume 
Percent * 



8,600 Original Res. Pressure 
8,260 Dew Point Pressure 

7.500 

6.500 

5.500 

4.500 
3,000 

1.500 



0 


0 


0 


0 


14.09 


2.48 


31.55 


5.55 


42.33 


7.45 


47.03 


8.27 


48.64 


8.56 


45.57 


6.68 



* Percent of Reservoir Hydrocarbon Pore Space 



HYDROCARBON ANALYSES OF GAS AND LIQUID REMAINING IN RESERVOIR 
AFTER DEPLETION TO 1,500 PSIA at 240°F 



Component 


Reservoir 


Gas 


Reservoir 


Liquid 


Mol % 15 


GPM@ 

.025 PSIA 


Mol 7. 


Liq . 


Vol . 7. 


Carbon Dioxide 


0.74 










Methane and Nitrogen 


92.36 




26.74 


7. 


42 


Ethane 


3. 70 




3. 22 


1 . 


34 


Propane 


1.46 


0.410 


1.74 


0 . 


79 


Iso-Butane 


0.44 


0.147 


0.65 


0 . 


35 


N-Butane 


0.31 


0.100 


0.65 


0 . 


34 


Iso-Pentane 


0.17 


0.063 


0.98 


0 . 


59 


N-Pentane 


0.11 


0.041 


0.74 


0 . 


44 


Hexanes 


0.21 


0.088 


1.99 


1 . 


34 


Heptanes Plus 


0.50 


0.251 


63.29 


87. 


39 


Totals 


100.00 


1.100 


100.00 


100. 


00 


Properties of Heptanes 


Plus : 










API Gravity = 


38.1° @ 60°F 










Specific Gravity = 


0.8343 










Molecular Weight = 


222 











DEPLETION STUDY AT 240°F 

HYDROCARBON ANALYSES OF PRODUCED WELL STREAM (MOL PER CENT) 
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Component 


Dew Point 
Fluid 

8260 




Reservoir Pressure - 


PSIA 




7500 


6500 


5500 


4500 


3000 


1500 


Carbon Dioxide 


0.67 


0.73 


0.84 


0.78 


0.75 


0.78 


0.74 


Methane and Nitrogen 


90.20 


91.08 


91.92 


92.53 


92.86 


92.88 


92.36 


Ethane 


3.56 


3.41 


3.29 


3.20 


3.24 


3.31 


3.70 


Propane 


1.37 


1.24 


1.12 


1.15 


1.15 


1.27 


1.46 


Iso-Butane 


0.43 


0.39 


0.36 


0.35 


0.35 


0.39 


0.44 


N-Butane 


0.35 


0.31 


0.27 


0.26 


0.26 


0.27 


0.31 


Iso-Pentane 


0.17 


0.17 


0.16 


0.16 


0.16 


0.15 


0.17 


N-Pentane 


0.13 


0.12 


0.12 


0.12 


0.11 


0.11 


0.11 


Hexanes 


0.30 


0.27 


0.23 


0.21 


0.20 


0.19 


0.21 


Heptanes Plus 


2.82 


2.28 


1.69 


1.24 


0.92 


0.65 


0.50 


Totals 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 



Properties of Heptanes Plus 

Specific Gravity 0.8256 

Molecular Weight 202 

GPM Content of Produced Well 


0.8198 

196 

Stream 


0.8146 0.8058 

186 169 

(Gal/MSCF) 


0.7981 

156 


0.7826 

137 


0.7665 

119 


Propane 


0.385 


0.349 


0.315 


0.323 


0.323 


0.357 


0.410 


Iso-Butane 


0.144 


0.130 


0.120 


0.117 


0.117 


0.130 


0.147 


N-Butane 


0.113 


0.100 


0.087 


0.084 


0.084 


0.087 


0.100 


Iso-Pentane 


0.063 


0.063 


0.060 


0.060 


0.060 


0.056 


0.063 


N-Pentane 


0.048 


0.044 


0.044 


0.044 


0.041 


0.041 


0.041 


Hexanes 


0.126 


0.113 


0.096 


0.088 


0.084 


0.080 


0.088 


Heptanes Plus 


2.233 


1.764 


1.249 


0.841 


0.582 


0.368 


0.251 


Totals 


3. 112 


2.563 


1.971 


1.557 


1.291 


1.119 


1.100 


Deviation Factor "Z" 


1.280 


1.210 


1.126 


1.054 


0.992 


0.927 


0.932 


Well Stream Produced 


Cumulative Per Cent 


0 


4.10 


11.07 


19.57 


29.84 


49.56 


73.51 
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Pseudo-critical Pressure - 662.30 PSIA 
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+ *<!* <j" <-0 lO 

CM O O O vD O O 

CJ CJ CJ ^ CJ r-H 2: 




o 00 vD <J CO CM 1—1 



Mol % 



Pressure (PSIA) 




Relative Volume (V/Vsat) 
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PRESSURE -VOLUME RELATIONS OF RESERVOIR FLUID @ 240°F 




Gallons per Thousand Cubic Feet 
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GASOLINE CONTENT OF HYDROCARBONS PRODUCED 



DURING PRESSURE DEPLETION @ 240°F 




0 2000 4000 6000 8000 



Pressure (PSIA) 




Deviation Factor n Z 
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DEVIATION FACTOR "Z" OF GAS PHASE 
DURING DEPLETION @ 240°F 



o Calculated 



• Experimental 




4000 



Pressure (PSIA) 



0 



2000 



6000 



8000 
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VOLUME OF WELL STREAM PRODUCED 
DURING DEPLETION @ 240°F 




Pressure (PSIA) 




Volume of Retrograde Liquid 
Percent of Hydrocarbon Pore Space 
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RETROGRADE CONDENSATION DURING 
DEPLETION @ 240 °F 




0 2000 4000 6000 8000 

Pressure (PSIA) 




APPENDIX III 

VOLUMETRIC SWEEP-OUT PREDICTION 
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POTENTIAL DISTRIBUTION PROGRAM 
PROGRAM POTMAP 
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PROGRAM POTMAP 

DIMENSION P(33, 38), RK (33,381,0(33,38) 





DO 40 


1=1.31 


40 


READ 


62, ( P ( I , J ) ,J = 1 ,12) 




DO 42 


1=1,31 


42 


READ 


62 , ( P ( I , J ) , J=13 , 24 ) 




DO 63 


1=1,31 


63 


READ 


62, ( P ( I , J ) , J = 2 5 , 36 ) 




DO 46 


1=1,31 


46 


READ 


41, ( Q ( I ,J) 




DO 48 


1=1.31 


48 


READ 


41, ( Q ( I , J ) . J = 1 9 , 36 ) 




DO 43 


1=1,31 


43 


READ 


49, ( RK ( I ,J) , J= 1 , 1 2 ) 




DO 45 


1=1.31 


45 


READ 


49, ( RK ( I , J ) . J=1 3 ,24 ) 




DO 61 


1=1,31 


61 


READ 


49, (RK( I , J) , J = 2 5 , 36 ) 




PRINT 


8 


8 


FORMAT ( 1H1 ) 




DO 1 


1=1.31 


1 


PRINT 


31 , ( P ( I , J ) » J = 1 . 18 ) 




PRINT 


8 




DO 2 


1=1.31 


2 


PRINT 


31, (PU»J) »J = 19»36) 




PRINT 


8 




DO 3 


1=1,31 


3 


PRINT 


13, (RK( I ,J).,J = 1 ,18 ) 




PRINT 


8 




DO 4 


1=1.31 


4 


PRINT 


13, (RKt I ,J) , J= 1 9 , 36 ) 


13 


FORMAT (1X,18F5.2//) 




LF = 


200 




ERROR 


= 1.0 




DO 10 


N= 1 » LF 




fail 


= 0.0 


12 


DO 29 


1=2,30 




DO 29 


J = 2 , 3 5 




DO 15 


K = 1 » 1 3 


15 


P ( 10, 


K)=P( 11, K) 




P( 10, 


14) = P( 11.15) 



P ( 9 * 14 ) = P ( 9 * 1 5 ) 

P ( 8 * 1 4 ) = P ( 8 * 1 5 ) 
DO 19 K= 15*19 
P ( 7 , K ) = P ( 8 * K ) 

P ( 7 » 20 ) = P ( 8 * 2 1 ) 
P ( 6 » 2 0 ) = P ( 6 * 2 1 ) 
P ( 5,20) = P(6 ,21 ) 
P( 5*22) =P( 6*22) 

P( 5,23) =P(6*23 ) 

P( 5,21) = P ( 6 , 2 1 ) 
P( 5,24) =P ( 6 , 2 4 ) 

P( 5,25) = P ( 6 , 26 ) 
P ( 4 , 2 5 ) = P ( 4 , 2 6 ) 
P( 3,25 ) = P( 3,26 ) 
P( 2,25) = P( 3,26 ) 



19 
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P ( 2 » 26 ) = P I 3 » 26 ) 

P ( 2 » 2 7 ) = P ( 3 * 2 7 ) 

P I 2 , 28 ) = P( 3,28 ) 

P( 2,29) = P( 3 ,30 ) 

P( 1,29) =P( 1,30) 

DO 18 K = 1,17 
18 PI 30, K) =P( 29, K) 

P ( 29 » 18 ) =P ( 29 , 17 ) 

PI 30, 18 ) = PI 29, 17) 

PI 28, 18 ) = PI 28, 17) 

PI 27,18 ) = P I 26 , 17 ) 

DO 16 K= 1 9 » 26 

16 PI 27, K) = PI 26, K) 

PI 27,27) = P I 26 , 26 ) 

P I 2 6 , 27 ) = P I 26 , 26 ) 

PI 25,27) = P I 2 5 , 26 ) 

PI 24,27) = P I 23 , 2 6 ) 

DO 17 K = 27 ,36 

17 P ( 29 , K ) = PI 23, K) 

CONST = 106000. /PI I,J) 

IF (PI I , J) ) 51,29,51 

51 IF (RK( I ,J) >52,29,52 

520PN= ( -CONST *Q( I ,J)+ I RK ( I+1,J)+RK( I,J))*P(I+1,J)+(RK( I-1,J)+RK 
1 I I , J) >*P I 1-1 , J ) + I RK( I , J+l ) +RK I I , J ) ) *P I I , J+l ) + (RKI I , J-l )+RK( I , J ) ) * 
2P(I,J-1)) / I RKI 1+1 , J) +RKI 1-1 , J)+RK( I ,J+1 )+RK( I ,J-1 ) +9.* 

3RK I I , J) ) 

IF (ERROR-ABSFIPN-PI I ,J) ) ) 26,27,27 

26 FAIL = 1.0 

27 PI I ,J) = PN 

29 CONTINUE 
PRINT 31, FAIL 

IF (FAIL) 33,33,10 
10 CONTINUE 
PRINT 8 

33 DO 30 1=1,31 

30 PRINT 31, (PI I,J) , J= 1 » 1 8 ) 

31 FORMAT I IX, 18F5.0// ) 

PRINT 8 

DO 32 1=1,31 

32 PRINT 31, I P I I , J ! , J=19 ,36 ) 

62 FORMAT (12F6.0) 

99 FORMAT I 12F6.2 ) 

91 FORMAT (18F9.0) 

99 FORMAT I 18F9.3 ) 

END 



END 



POTMAP 
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FUTURE PERFORMANCE PREDICTION 
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FUTURE PREDICTION PROGRAM 
PROGRAM PREDICT 
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104 





STOP 



CALCULATION OF GAS DEVIATION FACTOR 
Subroutine Z CALC (S,Z) 




CALCULATION OF DIMENSIONLESS WATER INFLUX 
Subroutine QFUNCT (DELT, QT) 
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WATER INFLUX CONDITION 
Subroutine INFLUX (DELTD, WD) 
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S=S+5 .0 



S=S-5 . 0 



>0 



< 0 



S=S+1.0 



S=S-1.0 \ 



> 0 



S=S+0. 1 [■*- 



< 0 



s=s-o.i ! • 



>0 




PRESSURE CALCULATION 
SUBROUTINE PRESS (P L POZ, P 2 > 



( DIMENSION ) 

W 

C COMMON"') 



s = p l 


- 5.0 






Calculate Z 



I 



P0ZT = f 

Li 

HZ 



E = P0ZT - P0Z 




t 
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Ill 



800 

101 

102 



200 

1000 



201 

1010 

C 



202 

1020 

C 



203 



204 



205 

1030 

C 



PROGRAM PREDICT 

DIMENSION A(6>6) , B ( 6 ) » T ( 47 1 ) , Q ( 47 1 ) ,DELP(999) ♦P(40) 

COMMON A,B»TEMP»PC»TC»T , Q » K , DEL D 
PRINT 800 
FORMAT ( 1 H 1 ) 

PRINT 101 

FORMAT ( 48H RESERVOIR COMPUTER PROJECT FUTURE PERFORMANCE,//) 

PRINT 102 

FORMAT ( 19H RESERVOIR DATA FOR,/) 

DATA INPUT FOR SUBROUTINE ZCALC CONSTANTS 
DO 1000 J= 1 ,6 , 1 
DO 1000 1=1, 6,1 
READ 200, A( I , J) 

FORMAT ( 13X, F14.ll ) 

CONTINUE 
DO 1010' 1 = 1 ,6 , 1 
READ 201 ,B( I ) 

FORMAT! F14. 11 ) 

CONTINUE 

DATA INPUT FOR SUBROUTINE QFUNCT TABLE 
DO 1020 1=1,471,3 

READ 202, T(I ) ,Q( I) ,T( I + 1),Q(I + 1) » T ( I +2 ) ,0(1+2) 

FORMAT! F12. 2, F 12. 3, FI 2.2, F 12. 3, FI 2. 2, FI 2.3 ) 

CONTINUE 

DATA INPUT FOR RESERVOIR PROPERTIES 

READ 203,AIWEL 

READ 203 , CONV 

READ 203, RI 

READ 203, DEPTH 

READ 203, TEMP 

READ 203, SGI 

READ 203,0 

READ 203, AH 

READ 203, SG 

READ 203, U 

READ 203, AK 

READ 203, TIME 

READ 203, SGR 

READ 203, SWC 

READ 203, PI 

READ 203, PC 

READ 203, TC 

READ 203, BW 

READ 203 , CPI 

READ 203, WIC 

READ 203, G 

FORMAT! F16. 7 ) 

DATA INPUT FOR PAST PERFORMANCE 
READ 204, GAS, DIST, WATER 
FORMAT! El 6.8, 4X, El 6.8 ,4X,E 16.8 ) 

I NC = CP I 

DO 1030 1 = 1 , INC, 1 
READ 205, P( I ) 

FORMAT! F10. 3 ) 

CONTINUE 

DATA INPUT FOR FUTURE PERFORMANCE 
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READ 700, KK 
FORMAT! 13 ) 

READ 203, GR 
READ 203, DR 
READ 203, WR 

DATA PRINT OUT FOR RESERVOIR PROPERTIES 
PRINT 300, INC 

FORMAT ( 45H INC, CURRENT NUMBER OF PRODUCTION I NCREMENTS= , 1 4 ) 

PRINT 3 01 , A I WEL 

FORMAT ( 48H MAX. NUMBER OF PRODUCING WELLS AT ANY ONE TIME=,F6.1) 
PRINT 302, BW 

FORMAT ( 50H BW, FORMATION VOLUME FACTOR FOR WATER IN BBLS /STB= , F6. 3 ) 
PRINT 303 , CONV 

FORMAT ( 44H CONV,SCF OF VAPOR PER BARREL OF D I ST I LLATE= , F7. 1 ) 

PRINT 304, RI 

FORMAT ( 48H RI, RADIUS TO INITIAL EDGE OF RESERVOIR IN FEET=,F8.1) 
PRINT 305, DEPTH 

FORMAT ( 34H DEPTH, DEPTH OF RESERVOIR IN FEET=,F8.1) 

PRINT 306, TEMP 

FORMAT ( 50H T EMP , TEMPER A TURE OF RESERVOIR IN DEGREES RANK I NE= , F7 . 1 ) 
PRINT 307, PC 

FORMAT ( 36H PC , PSEUDOCR I T I CAL PRESSURE IN PSIA=,F6.1) 

PRINT 308, TC 

FORMAT ( 50H T C , PS EUDOC R I T I C AL TEMPERATURE IN DEGREES RANK I NE= , F6 . 1 ) 
PRINT 309 , SGI 

FORMAT ( 40H SGI, INITIAL GAS SATURATION AS FR AC T I ON= , F 5 . 3 ) 

PRINT 310,0 

FORMAT ( 24H 0, POROSITY AS FR ACT I ON= , F 5 . 3 ) 

PRINT 311, AH 

FORMAT ( 43H AH, AVERAGE THICKNESS OF RESERVOIR IN FEET=,F6.1) 

PRINT 312, SG 

FORMAT ( 44H SG, SPECIFIC GRAVITY OF GAS RELATIVE TO AIR=,F7.4> 

PRINT 3 1 3 , U 

FORMAT ( 46H U, VISCOSITY OF RESERVOIR WATER IN CENT I PO I SE= , F6 • 3 ) 
PRINT 314, AK 

FORMAT ( 47H AK, WATER PERMEABILITY OF RESERVOIR IN DARC I ES = , F7 . 4 ) 
PRINT 315, TIME 

FORMAT ( 33H TIME, TIME INTERVAL USED IN DAYS=,F7.2) 

PRINT 316, SGR 

FORMAT ( 4 1H SGR, RESIDUAL GAS SATURATION AS FR ACT I ON = , F5 . 3 ) 

PRINT 317, SWC 

FORMAT ( 42H SWC, CONNATE WATER SATURATION AS FRACT I ON= , F5 * 3 ) 

PRINT 318, WIC 

FORMAT ( 45H WIC, WATER INFLUX CONSTANT IN BBLS. PER PS I A= , F 1 0 . 3 , / / ) 
DATA PRINT OUT FOR RESERVOIR PRODUCTION 
PRINT 103 

FORMAT ( 55H RESERVOIR PRODUCTION FOR FUTURE PERFORMANCE PREDICTION) 
PRINT 104 
FORMAT ( 1H ) 

PRINT 319, GR 

FORMAT ( 41H AVERAGE DAILY DRY GAS PRODUCTION IN SCF=,F16.1) 

PRINT 320, DR 

FORMAT ( 44H AVERAGE DAILY DISTILLATE PRODUCTION IN STB=*F16.1) 

PRINT 321, WR 

FORMAT ( 39H AVERAGE DAILY WATER PRODUCTION IN STB= , F16. 1 , // ) 

MAIN BODY OF PROGRAM 
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105 

322 

323 

400 

401 

106 
107 

2 



3 



402 




S=P( 1 ) 

CALL ZCALC ( S > Z ) 

C=0.00514590*TEMP 

BGI=C*Z/S 

DELTD=1858064.0*AK*T I ME/ (U*0*RI*RI ) 

GRB=G*BGI 

FF=( 1.0-SGR-SWC) /SGI 

FLOOD=GRB*FF 

PRINT 105 

FORMAT ( 1 9H INITIAL CONDITIONS,/) 

PRINT 322 , PI 

FORMAT ( 26H INITIAL PRESSURE IN PSIA=,F8.2) 

PRINT 323, G 

FORMAT ( 29H INITIAL GAS IN PLACE IN SCF=»E16.8) 

PRINT 400, FLOOD 

FORMAT ( 47H MAXIMUM FLOODABLE RESERVOIR VOLUME IN BARRELS= , E16 . 8 ) 
PRINT 401, Z 

FORMAT ( 36H INITIAL GAS COMPRESSIBILITY F ACTOR= , E 1 6 . 8 , / / ) 

PRINT 106 

FORMAT ( 1 H »/,18H FUTURE PREDICTION,/) 

PRINT 107 

FORMAT { 1 5H TIME PRESSURE, 9X,2HWE,16X, 3 HGPW,8X,7HDELTA P) 

DELPt 1)=0.0 
DO 2 1=2, INC, 1 
DELPf I )=P( 1-1 ) -P ( I ) 

P1=P( INC) 

WT=WATER 
DT = D I ST 
GT = GAS 
W= WR*T I ME 
D=DR*TIME*CONV 
GDP=GR*TIME 
K= INC 
K = K + 1 
L = K- 1 



WT = W+WT 

DT = D+DT 

GT =GDP+GT 

EKAY=GT +DT 

DELP( K) =0.0 

CALL INFLUX(DELTD,WD) 

WE=W I C*WD 

POZ=C*( G-EKAY ) / ( GR8+WT*BW-WE ) 

CALL PRESS(P1,P0Z,P2) 

DELP(K) =P1-P2 

PRINT 402,L,P2,WE,EKAY,DELP(K) 

FORMAT ( I5,2X,F8.2,2X,E16.8,2X,E16.8,2X,F8.2) 

P1 = P2 

I F ( WE-FLOOD ) 8,4,4 
PRINT 108 

FORMAT ( 1H » / , 2 7H ALL WELLS HAVE FLOODED OUT,/////) 

CONTINUE 

IF(KK-l) 6,6,7 

CONTINUE 

GO TO 1 

IF ( P 1-14 • 7 ) 10,10,9 
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9 GO TO 3 
10 PRINT 109 

109 FORMAT ( 1H */.19H RESERVOIR DEPLETED,/////) 
CONTINUE 
GO TO 5 
END 
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SUBROUTINE ZCA|_C(S,Z) 

DIMENSION A(6»6)»B(6)»T( 471) >0(471) >DELP(999) >P(40) 
DIMENSION PX ( 6 ) »PY(6) 

COMMON A»B»TEMP»PC»TC>T »Q > K. » DELP 
Y= ( ( 2.0*TEMP/TC) -4.0 ) /I .9 
P Y ( 1 ) =B ( 1 ) 

PY ( 2 ) =B ( 2 )*Y 

PY(3)=B(3)*( (3.0*Y*Y)-1.0) 

PY (4 ) =B ( 4 )* I ( 5 » 0*Y*Y ) ~3 ♦ 0 ) *Y 
PY(5)=B(5)*(( ( 35 .0*Y*Y )-30.0)*Y*Y+3.0) 

PY ( 6 ) =B ( 6 )*( ( I (63.0*Y*Y)-70.0)*Y*Y)+15.0)*Y 
X=( (2. 0*5 /PC 1-15.0) /14.8 
PX ( 1 ) =B I 1 ) 

PX ( 2 ) =B ( 2 )*X 

PX ( 3 ) =B ( 3 ) * ( ( 3 • 0*X*X ) - 1 • 0 ) 

PX(4)=B(4)*( (5.0*X*X)-3.0)*X 

PX ( 5 ) =8 ( 5 ) * ( ( ( ( 35 • 0*X*X )-30.0)*X*X)+3.0) 

PX(6)=B(6)*(((( 63 . 0*X*X ) -70 . 0 ) *X*X ) +15 . 0 ) *X 

APP = 0 . 0 

DO 1 M=1 »6 »1 

DO 1 N= 1 » 6 » 1 

1 APP=( A(N,M)*PX(M)*PY(N) )+APP 
Z= APP 
RETURN 
END 
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SUBROUTINE QFUNCT (DELT »QT ) 

DIMENSION A ( 6 » 6 ) » B ( 6 ) » T ( 47 1 ) » Q ( 47 1 ) * DELP ( 999 ) , P(40 ) 
COMMON A*B» TEMP, PC,TC,T ,Q,K, DELP 
TT = DELT 

DO 3 1=1*471, 1 
I F ( TT-T ( I ) ) 1,2,3 

1 DQ 1= ( Q ( I - 1 ) -Q (1-2) )/(T(I-l)-T{ 1-2) ) 

DQ2 = ( ( Q ( I ) — Q { I - 1 ) )/(T( I ) — T ( I — 1 ) )-DQl ) / ( T ( I ) -T ( 1-2 ) ) 
QT=Q( 1-2 )+( TT-T ( 1-2) )*DQ1+(TT-T( 1-2) )*( TT-T ( 1-1 ) )*DQ2 
GO TO 4 

2 QT=Q( I ) 

GO TO 4 

3 CONTINUE 

4 RETURN 
END 
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SUBROUTINE I N FLUX ( DEL T D , WD ) 

DIMENSION A ( 6 > 6 ) > B ( 6 ) » T ( 47 1 ) » Q ( A7 1 ) >DELP(999) >P(40) 
COMMON A>B»TEMP»PC»TC»T »Q?K»DELP 
WD=0 . 0 

DO 1 J= 1 » K ♦ 1 
T I = J- 1 

DELT=TI*DELTD 

CALL QFUNCT ( DELT »QT ) 

N=K+1-J 
W I =DELP ( N ) *QT 
1 WD=W I +WD 
CONTINUE 
RETURN 
END 
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SUBROUTINE PR ESS ( P 1 , POZ , P2 ) 

DIMENSION A ( 6 * 6 ) » B ( 6 ) » T ( 47 1 ) » 0 ( 47 1 ) »DELP(999) ,P(40) 
COMMON A,B,TEMP,PC,TC,T ,Q,K,DELP 
S= P 1-5 • 0 
CALL ZCALC(S,Z) 

POZT=S/Z 

E=POZT-POZ 

IF(ABSF(E)-5.0) 5*5,2 

I F ( E ) 4,17,3 

S=S-5 • 0 

GO TO 1 

S = S + 5 • 0 

GO TO 1 

I F ( ABSF ( E ) -1 • 0 ) 9,9,6 

I F ( E ) 8,17,7 

S=S-1 .0 

GO TO ' 1 

S=S+1 .0 

GO TO 1 

I F ( ABSF ( E ) -0 « 1 ) 13,13,10 

I F ( E ) 12,17,11 

S=S-0.1 

GO TO 1 

S=S+0.1 

GO TO 1 

IF(ABSF(E)-0.02) 17,17,14 

IF ( E) 16,17,15 

S=S-0.01 

GO TO 1 

S=S+0.01 

GO TO 1 

P2 = S 

RETURN 

END 



END 



PREDICT 
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C DATA INPUT FOR SUBROUTINE ZCALC CONSTANTS 



2 • 1433504 
0.083176184 
-0.021467042 
-0.00087140318 
0.0042846283 
-0.0016595343 
0.33123524 
-0.13403614 
0.066880961 
-0.027174261 
0.0088512291 
-0.0021520929 
0.10572871 
-0.050393654 
0.0050924798 
0.010551336 
-0.0073181933 
0.0026959963 
-0.052184040 
0.044312146 
-0.019329465 
0.0058972516 
0.0015366676 
-0.0028326809 
0.019703980 
-0.026383354 
0.019262143 
-0.011535390 
0.0042910089 
-0.00081302526 
-0.0053095900 
0.0089178330 
-0.010894821 
0.009559389 
-0.0060114017 
0.0031175170 

0.7071068 

1.224745 

0.7905695 

0.9354145 

0.265165 

0.293151 



DATA INPUT 


FOR SUBROUTINE QFUNCT 


TABLE 






00000000 


00000000 


oooooooi 


00000112 


00000005 


00000278 


00000010 


00000404 


00000015 


00000520 


00000020 


00000606 


00000025 


00000689 


00000030 


00000758 


00000040 


00000898 


00000050 


00001020 


00000060 


00001140 


00000070 


00001251 


00000080 


00001359 


00000090 


00001469 


00000100 


00001569 


00000200 


00002447 


00000300 


00003202 


00000400 


00003893 



120 



00000500 

00000800 

00001100 

00001400 

00001700 

00002000 

00002300 

00002600 

00002900 

00003200 

00003500 

00003800 

00004100 

00004400 

00004700 

00005000 

00005300 

00005600 

00005900 

00006200 

00006500 

00006800 

00007100 

00007400 

00007700 

00008000 

00008300 

00008600 

00008900 

00009200 

00009500 

00009800 

00010500 

00012000 

00013500 

00015000 

00016500 

00018000 

00019500 

00021000 

00022500 

00024000 

00025500 

00027000 

00028500 

00030000 

00031500 

00033000 

00034500 

00036000 

00037500 

00039000 

00040500 

00042000 

00043500 

00045000 



00004539 
00006314 
00007940 
00009461 
00010913 
00012319 
00013684 
00015013 
00016313 
00017590 
00018845 
00020080 
00021298 
00022500 
00023684 
00024855 
00026020 
00027174 
00028314 
00029443 
00030565 
00031679 
00032785 
00033883 
00034974 
00036058 
00037136 
00038207 
00039272 
00040331 
00041385 
00042433 
00044858 
00049968 
00054976 
00059985 
00064737 
00069512 
'"'00 742 2 6 
00078886 
00083497 
00088062 
00092589 
00097081 
00101540 
00105968 
00110367 
00114738 
00119083 
00123403 
00127699 
00131972 
00136223 
00140453 
00144664 
00148856 



00000600 
00000900 
00001200 
00001500 
00001800 
00002100 
00002400 
00002700 
00003000 
00003300 
00003600 
00003900 
0 0004200 
00004500 
00004800 
00005100 
00005400 
00005700 
00006000 
00006300 
00006600 
00006900 
00007200 
00007500 
00007800 
00008100 
00008400 
00008700 
00009000 
00009300 
00009600 
00009900 
00011000 
00012500 
00014000 
00015500 
00017000 
00018500 
00020000 
00021500 
00023000 
00024500 
00026000 
00027500 
00029000 
00030500 
00032000 
00033500 
00035000 
00036500 
00038000 
00039500 
00041000 
00042500 
00044000 
00045500 



00005153 

00006869 

00008457 

00009949 

00011386 

00012778 

00014131 

00015450 

00016742 

00018011 

00019259 

00020488 

00021701 

00022897 

00024076 

00025244 

00026406 

00027555 

00028691 

00029818 

00030937 

00032048 

00033151 

00034247 

00035336 

00036418 

00037494 

00038563 

00039626 

00040684 

00041735 

00042781 

00046574 

00051648 

00056625 

00061517 

00066336 

00071090 

00075785 

00080428 

00085023 

00089575 

00094090 

00098571 

00103019 

00107437 

00111827 

00116189 

00120526 

00124838 

00129126 

00133391 

00137635 

00141859 

00146064 

00150249 



00000700 

00001000 

00001300 

00001600 

00001900 

00002200 

00002500 

00002800 

00003100 

00003400 

00003700 

00004000 

00004300 

00004600 

00004900 

00005200 

00005500 

00005800 

00006100 

00006400 

00006700 

00007000 

00007300 

00007600 

00007900 

00008200 

00008500 

00008800 

00009100 

00009400 

00009700 

00010000 

00011500 

00013000 

00014500 

00016000 

00017500 

00019000 

00020500 

00022000 

00023500 

00025000 

00026500 

00028000 

00029500 

00031000 

00032500 

00034000 

00035500 

00037000 

00038500 

00040000 

00041500 

00043000 

00044500 

00046000 



00005743 

00007411 

00008964 

00010434 

00011855 

00013233 

00014573 

00015883 

00017167 

00018429 

00019671 

00020894 

00022101 

00023291 

00024466 

00025633 

00026791 

00027935 

00029068 

00030192 

00031308 

00032417 

00033517 

00034611 

00035697 

00036777 

00037851 

00038919 

00039979 

00041034 

00042084 

00043129 

00048277 

00053317 

00058265 

00063131 

00067928 

00072661 

00077338 

00081965 

00086545 

00091084 

00095588 

00100057 

00104495 

00108904 

00113284 

00117638 

00121966 

00126270 

00130550 

00134808 

00139045 

00143262 

00147461 

00151640 
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00046500 

00048000 

00049500 

00052000 

00054000 

00057000 

00059000 

00062000 

00064000 

00067000 

00069000 

00072000 

00074000 

00077000 

00079000 

00082000 

00084000 

00087000 

00089000 

00092000 

00094000 

00097000 

00099000 

00102000 

00104000 

00107000 

00109000 

00112000 

00114000 

00117000 

00119000 

00122000 

00124000 

00127000 

00129000 

00132000 

00134000 

00137000 

00139000 

00142000 

00144000 

00147000 

00149000 

00155000 

00162500 

00170000 

00177500 

00185000 

00192500 

00200000 

00207500 

00215000 

00222500 

00230000 

00237500 

00245000 



00153029 

00157184 

00161322 

00168183 

00173639 

00181774 

00187166 

00195208 

00200542 

00208502 

00213784 

00221670 

00226904 

00234721 

00239912 

00247668 

00252819 

00260515 

00265629 

00273274 

00278353 

00285498 

00290995 

00298543 

00303560 

00311066 

00316055 

00323517 

00328480 

00335906 

00340843 

00348230 

00353144 

00360496 

00365386 

00372704 

00377572 

00384859 

00389705 

00396959 

00401786 

00409013 

00413820 

00428196 

00446077 

00463863 

00481558 

00499167 

00516695 

00534145 

00551522 

00568830 

00586072 

00603252 

00620372 

00637437 



00047000 

00048500 

00050000 

00052500 

00055000 

00057500 

00060000 

00062500 

00065000 

00067500 

00070000 

00072500 

00075000 

00077500 

00080000 

00082500 

00085000 

00087500 

00090000 

00092500 

00095000 

00097500 

00100000 

00102500 

00105000 

00107500 

''>0110000 

00112500 

00115000 

00117500 

00120000 

00122500 

00125000 

00127500 

00130000 

00132500 

00135000 

00137500 

00140000 

00142500 

00145000 

00147500 

00150000 

00157500 

00165000 

00172500 

00180000 

00187500 

00195000 

00202500 

00210000 

00217500 

00225000 

00232500 

00240000 

00247500 



00154416 

00158565 

00162698 

00169549 

00176357 

00183124 

00189852 

00196544 

00203201 

00209825 

00216417 

00222980 

00229514 

00236020 

00242501 

00248957 

00255388 

00261795 

00268181 

00274545 

00280888 

00287211 

00293514 

00299799 

00306065 

00312314 

00318545 

00324760 

00330958 

00337142 

00343308 

00349460 

00355597 

00361720 

00367828 

00373922 

00380003 

00386070 

00392125 

00398167 

00404197 

00410214 

00416220 

00434168 

00452016 

00469771 

00487437 

00505019 

00522520 

00539945 

00557299 

00574585 

00591806 

00608965 

00626066 

00643113 



0 0047500 
00049000 
00051000 
00053000 
00056000 
00058000 
00061000 
00063000 
00066000 
00068000 
00071000 
00073000 
00076000 
00078000 
00081000 
00083000 
00086000 
00088000 
00091000 
00093000 
00096000 
00098000 
00101000 
00103000 
00106000 
00108000 
00111000 
00113000 
00116000 
00118000 
00121000 
00123000 
00126000 
00128000 
00131000 
00133000 
00136000 
00138000 
00141000 
00143000 
00146000 
00148000 
00152500 
00160000 
00167500 
00175000 
00182500 
00190000 
00197500 
00205000 
00212500 
00220000 
00227500 
00235000 
00242500 
00250000 



00155801 

00159945 

00165444 

00170914 

00179069 

00184473 

00192533 

00197878 

00205854 

00211145 

00219046 

00224289 

00232120 

00237218 

00245086 

00250245 

00257953 

00263073 

00270729 

00275815 

00283420 

00288473 

00296030 

00301053 

00308567 

00313562 

00321032 

00326000 

00333433 

00338376 

00345770 

00350688 

00358048 

00362942 

00370267 

00375139 

00382432 

00387283 

00394543 

00399373 

00406606 

00411418 

00422214 

00440128 

00457945 

00475669 

00493307 

00510861 

00528337 

00545737 

00563068 

00580332 

00597532 

00614672 

00631755 

00648781 
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00255000 

00270000 

00285000 

00300000 

00315000 

00330000 

00345000 

00360000 

00375000 

00390000 

00405000 

00420000 

00435000 

00450000 

00465000 

00480000 

00495000 

00520000 

00550000 

00580000 

00610000 

00640000 

00670000 

00700000 

00730000 

00760000 

00790000 

00820000 

00850000 

00880000 

00910000 

00940000 

00970000 

01000000 

01750000 

03000000 

05000000 

07500000 

10000000 



00660093 

00693877 

00727449 

00760833 

00794042 

00827088 

00859974 

00892712 

00925309 

00957773 

00990108 

01022318 

01054409 

01086390 

01118257 

01150012 

01181666 

01234203 

01296893 

01359225 

01421224 

01482912 

01544305 

01605418 

01666265 

01726859 

01787212 

01847336 

01907243 

01966943 

02026438 

02085744 

02144878 

02203861 

03633368 

05899508 

09363099 

13531457 

17586284 



00260000 

00275000 

00290000 

00305000 

00320000 

00335000 

003500O0 

00365000 

00380000 

00395000 

00410000 

00425000 

00440000 

00455000 

00470000 

00485000 

00500000 

00530000 

00560000 

00590000 

00620000 

00650000 

00680000 

00710000 

00740000 

00770000 

00800000 

00830000 

00860000 

00890000 

00920000 

00950000 

00980000 

01250000 

02000000 

03500000 

06000000 

08000000 

12500000 



00671379 

00705090 

00738598 

00771922 

00805075 

00838067 

00870903 

00903594 

00936144 

00968566 

01000858 

01033028 

01065082 

01097024 

01128854 

01160574 

01192198 

01255141 

01317709 

01379927 

01441820 

01503408 

01564706 

01625729 

01686490 

01747002 

01807278 

01867329 

01927166 

01986796 

02046227 

02105473 

02164555 

02688967 

04095800 

06780247 

11047299 

14350121 

21560732 



00265000 

00280000 

00295000 

00310000 

00325000 

00340000 

00355000 

00370000 

00385000 

00400000 

00415000 

00430000 

00445000 

00460000 

00475000 

00490000 

00510000 

00540000 

00570000 

00600000 

00630000 

00660000 

00690000 

00720000 

00750000 

00780000 

00810000 

00840000 

00870000 

00900000 

00930000 

00960000 

00990000 

01500000 

02500000 

04000000 

07000000 

09000000 

15000000 



00682640 

00716280 

00749725 

00782992 

00816090 

00849028 

00881816 

00914459 

00946966 

00979344 

01011595 

01043724 

01075743 

01107646 

01139439 

01171125 

01213222 

01276037 

01338486 

01400593 

01462383 

01523872 

01585077 

01646011 

01706688 

01767120 

01827319 

01887298 

01947065 

02006628 

02065996 

02125184 

02184216 

03164780 

05005726 

07650096 

12708358 

15975389 

25380000 



C 

2.0 
768.0 
2190. 
11450. 
690. 
.80 
.28 
64.4 
.7975 
.30 
. 050 
180 . 
.30 



DATA INPUT FOR RESERVOIR PROPERTIES 
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0.15 

8838. 

662.26 

387.48 

1.03 

7.0 

7.52757 

606898622000. 

2.07907540E+10 1 . 542 2 1 800E+06 

8838. 

8764. 

8657. 

8509. 

8364. 

8219. 

8100. 

9 

25000000. 

1250. 

30. 



8.67600000F+03 
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RESERVOIR COMPUTER PROJECT FUTURE PERFORMANCE 



RESERVOIR DATA FOR 

INOCURRENT NUMBER OF PRODUCTION I NCR EMENTS= 7 

MAX. NUMBER OF PRODUCING WELLS AT ANY ONE T I ME= 2.0 

BW* FORMATION VOLUME FACTOR FOR WATER IN BBLS/STB= 1.030 

CONV.SCF OF VAPOR PER BARREL OF D I ST I LLATE= 768.0 

RI. RADIUS TO INITIAL EDGE OF RESERVOIR IN FEET= 2190.0 

DEPTH , DEPTH OF RESERVOIR IN FEE T = 11450.0 

TEMP, TEMPERATURE OF RESERVOIR IN DEGREES RANK.I NE= 690.0 

PC, PSEUDO CRITICAL PRESSURE IN PSIA= 662.3 

TC,PSEUDOCRI TICAL TEMPERATURE IN DEGREES RANK I NE= 387.5 

SGI, INITIAL GAS SATURATION AS FRACTION= .800 

0, POROSITY AS FRACT I ON= .280 

AH, AVERAGE THICKNESS OF RESERVOIR IN FEET= 64.4 

SG, SPECIFIC GRAVITY OF GAS RELATIVE TO AIR= .7975 

U, VISCOSITY OF RESERVOIR WATER IN CENTIPOISE= .300 

AK, WATER PERMEABILITY OF RESFRVOIR IN DARC TES= .0500 

TIME, TIME INTERVAL USED IN DA YS= 180.00 

SGR, RESIDUAL GAS SATURATION AS FRACT I ON= .300 

SWC, CONNATE WATER SATURATION AS FRACTION^ .150 

WIC, WATER INFLUX CONSTANT IN'BBLS. PER PSIA= 7.528 



RESERVOIR PRODUCTION FOR FUTURE PERFORMANCE 
AVERAGE DAILY DRY GAS PRODUCTION IN SCF= 
AVERAGE DAILY DISTILLATE PRODUCTION IN STB= 
AVERAGE DAILY WATER PRODUCTION IN STB= 



PREDICTION 

25000000.0 

1250.0 

30.0 



INITIAL CONDITIONS 

INITIAL PRESSURE IN PSIA= 8838.00 

INITIAL GAS IN PLACE IN SCF= 6.06898622E 11 

MAXIMUM FLOODABLE RESERVOIR VOLUME IN BARRELS= 2.22914288E 08 
INITIAL GAS COMPRESSIBILITY FACTOR= 1.32981926E 00 



FUTURE PREDICTION 



TIME 


PRESSURE 


WE 




GPW 




DELTA P 


7 


8009.66 


3.01997662E 


05 


2 . 54650962E 


10 


90.34 


8 


7870.52 


3 . 92764487E 


05 


3.01378962E 


10 


139.14 


9 


7734.90 


4.9871 1813E 


05 


3.48 106962E 


10 


135.62 


10 


7602.48 


6. 172 96 089E 


05 


3 . 94834962E 


10 


132.42 


11 


7473.08 


7.47733989E 


05 


4.41 562962E 


10 


129.40 


12 


7346.57 


8.89432311E 


05 


4.88 290962E 


10 


126.51 


13 


7222.84 


1 .04188907E 


06 


5.35018962E 


10 


123.73 


14 


7101.81 


1 .20466079E 


06 


5.81 746962E 


10 


121.03 


15 


6983.39 


1.37735030E 


06 


6 • 28474962E 


10 


118.42 


16 


6867. 53 


1 . 55959584E 


06 


6 • 75202962E 


10 


115.86 


17 


6754. 18 


1 .75105577E 


06 


7. 21930962E 


10 


113.35 



18 

19 

20 

21 

22 

23 

24 

23 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 
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6643.27 


1.95 14067 6 E 


06 


7.68658962E 


10 


110.91 


6534.76 


2. 160348 59E 


06 


8.15386962E 


10 


108.51 


6428.60 


2.37757163E 


06 


8.621 14962E 


10 


106.16 


6324.75 


2.60283387E 


06 


9.08842962E 


10 


103.85 


6223. 15 


2.83585730E 


06 


9 . 55570962E 


10 


101.60 


6123.75 


3.0 76406 7 5 E 


06 


1 .00229896E 


11 


99.40 


6026.50 


3. 32422979E 


06 


1 .04902696E 


11 


97.25 


5931.35 


3.57910580E 


06 


1 .09575496E 


11 


95.15 


5838.24 


3.84080861E 


06 


1 . 14248296E 


11 


93.11 


5747.13 


4.10910827E 


06 


1 . 18921096E 


11 


91.11 


5657.94 


4.38383777E 


06 


1 .23593896E 


11 


89.19 


5570.63 


4.66479429E 


06 


1 .28266696E 


11 


87.31 


5485.14 


4.95179404E 


06 


1 . 32939496E 


11 


85.49 


5401.40 


5.24466308E 


06 


1 .37612296E 


11 


83.74 


5319.37 


5 . 54323484E 


06 


1 .42285096E 


1 1 


82.03 


5238.99 


5.84734876E 


06 


1 .46957896E 


11 


80.38 


5160. 19 


6. 15685365E 


06 


1 .51630696E 


11 


78.80 


5082.93 


6.47160388E 


06 


1 .56303496E 


11 


77.26 


5007. 15 


6.79146218E 


06 


1 .60976296E 


11 


75.78 


4932.79 


7.1 1 6 296 3 2 E 


06 


1 .65649096E 


11 


74.36 


4859.81 


7.44598239E 


06 


1 .70321896E 


11 


72.98 


4788.15 


7.78039908E 


06 


1 . 74994696E 


11 


71.66 


4717.76 


8.1 1943332E 


06 


1 . 79667496E 


11 


70.39 


4648 .59 


8.46297811E 


06 


1 .84340296E 


11 


69.17 


4580.60 


8. 81092942 E 


06 


1 .89013096E 


11 


67.99 


4513.74 


9. 16318972E 


06 


1 .93685896E 


11 


66.86 


4447.97 


9.51966598E 


06 


1 .98358696E 


11 


65.77 


4383.24 


9.88026859E 


06 


2.03031496E 


11 


64.73 


4319.51 


1.02449144E 


07 


2 . 07 7042 96E 


11 


63.73 


4256.74 


1 .06135227E 


07 


2.12 377096E 


11 


62.77 


4194.90 


1 .09860180E 


07 


2 . 17049896E 


11 


61.84 


4133.94 


1.13623269E 


07 


2.21722696E 


11 


60.96 


4073.83 


1 . 17423812E 


07 


2 . 26395496E 


11 


60.11 


4014.54 


1.21261164E 


07 


2.31068296E 


11 


59.29 


3956.03 


1.2513 46 99E 


07 


2 . 35741096E 


11 


58.51 


3898.27 


1.29043838E 


07 


2.40413896E 


11 


57.76 


3841.23 


1. 32988022E 


07 


2.45086696E 


11 


57.04 


3784.88 


1 . 36966727 E 


07 


2. 4975949 6E 


11 


56.35 


3729.18 


1.40979449E 


07 


2 . 54432296E 


11 


55.70 


3674. 12 


1.45025730E 


07 


2.59105096E 


11 


55.06 


3619.67 


1 • 49 1 05 1 1 2 E 


07 


2.63777896E 


11 


54.45 


3565.79 


1 .53217166E 


07 


2.68450696E 


11 


53.88 


3512.46 


1.57361517E 


07 


2.73 123496E 


11 


53.33 


3459.66 


1 . 6 1 53778 3 E 


07 


2.77796296E 


11 


52.80 


3407.38 


1.65745573E 


07 


2 . 82469096E 


11 


52.28 


3355.57 


1 .69984484E 


07 


2.87141896E 


11 


51.81 


3304.23 


1 .74254188E 


07 


2.91814696E 


11 


51.34 


3253.32 


1.7 8 5543 43 E 


07 


2 . 9648 7496E 


11 


50.91 


3202.82 


1 .82884669E 


07 


3.01 160296E 


11 


50.50 


3152.74 


1.8 724491 2 E 


07 


3.05833096E 


11 


50.08 


3103.03 


1 .91634781E 


07 


3. 10505896E 


11 


49.71 


3053.68 


1.9605404 IE 


07 


3.15 178696E 


11 


49.35 


3004.67 


2.00502457E 


07 


3. 19851496E 


11 


49.01 


2955.97 


2 . 049798 27E 


07 


3. 24524296E 


11 


48.70 


2907.60 


2.09485976E 


07 


3.29 197096E 


11 


48.37 


2859.50 


2. 14020682E 


07 


3 . 33869896E 


11 


48.10 



74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 
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109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 



126 



2811.68 2 

2764.11 2 

2716.78 2 

2669.68 2 

2622.78 2 

2576.07 2 

2529.54 2 

2483.18 2 

2436.96 2 

2390.86 2 

2344.88 2 

2299.05 2 

2253.25 2 

2207.58 2 

2161.95 2 

2116.35 2 

2070.82 2 

2025.29 3 

1979.77 3 

1934.24 3 

1888.69 3 

1843.09 3 

1797.46 3 

1751.76 3 

1705.96 3 

1660.09 3 

1614.10 3 

1568.00 3 

1521.77 3 

1475.36 3 

1428.76 3 

1382.01 3 

1335.01 3 

1287.81 3 

1240.39 3 

1192.68 3 

1144.68 4 

1096.38 4 

’1047.78 4 

* 998.79 4 

949.48 4 

899.74 4 

849.62 4 

799.02 4 

747.99 4 

' 696.45 4 

644.35 4 

591.74 4 

538.51 4 

484.66 4 

430.14 4 

374.92 4 

318.93 4 

262.20 5 

204.60 5 

146.12 5 



18583818E 


07 


23175214E 


07 


27794740E 


07 


3244226 5E 


07 


37117672E 


07 


41820867E 


07 


46551770E 


07 


5 1310297E 


07 


56096366E 


07 


60909930E 


07 


65750960E 


07 


70619420E 


07 


755152 16E 


07 


8 043845 5 E 


07 


85389026E 


07 


90366994E 


07 


95372392E 


07 


00405 1 93 E 


07 


0 546548 2 E 


07 


10553302E 


07 


15668727E 


07 


20811827E 


07 


25982721E 


07 


31181457E 


07 


36408157E 


07 


41662966E 


07 


4694596 5E 


07 


52257331E 


07 


57597186E 


07 


62965691E 


07 


68363065E 


07 


73789519E 


07 


79245192E 


07 


84730397E 


07 


90245323E 


07 


95790186E 


07 


01365305E 


07 


06970957E 


07 


12607427E 


07 


1 8275004E 


07 


23974101E 


07 


2970497 8E 


07 


35468084E 


07 


41263736E 


07 


47092427E 


07 


52954528E 


07 


58850553 E 


07 


64781051E 


07 


70746455E 


07 


7 674737 1 E 


07 


82784344E 


07 


88857997E 


07 


94968966 E 


07 


01117984E 


07 


07305695 E 


07 


1 3 5 3294 8 E 


07 



3. 38542696E 
3.43215496E 
3.47888296E 
3.52561096E 
3. 57233896E 
3.61 906696E 
3 . 66579496E 
3 • 71252296E 
3.75925096E 
3 . 80597896E 
3 . 85270696E 
3 • 89943496E 
3 .94616296E 
3 . 99289096E 
4. 03961896E 
4.08634696E 
4.13307496E 
4. 17980296E 
4. 22653096E 
4. 273258 96E 
4.31998696E 
4.36671496E 
4.4 1 3442 96E 
4.46017096E 
4. 50689896E 
4.55362696E 
4.60035496E 
4. 64708296E 
4.69381096E 
4.74053896E 
4.78726696E 
4.83399496E 
4.88072296E 
4.92745096E 
4. 97417896E 
5 . 02090696E 
5 . 06763496E 
5. 11436296E 
5.161 09096E 
5.20781896E 
5 • 2 5454696E 
5.30 127496E 
5 . 34800296E 
5.39473096E 
5.44145896E 
5.48818696E 
5.5349 1496E 
5.58 164296E 
5.62837096E 
5.67 509896E 
5.72 182696E 
5 . 76855496E 
5.81528296E 
5.86201096E 
5 . 90873896E 
5.95546696E 



11 


47.82 


11 


47.57 


11 


47.33 


11 


47.10 


11 


46.90 


11 


46.71 


11 


46.53 


11 


46.36 


11 


46.22 


11 


46.10 


11 


45.98 


11 


45.83 


11 


45.80 


11 


45.67 


11 


45.63 


11 


45.60 


11 


45.53 


11 


45.53 


11 


45.52 


11 


45.53 


11 


45.55 


11 


45.60 


11 


45.63 


11 


45.70 


11 


45.80 


11 


45.87 


11 


45.99 


11 


46.10 


11 


46.23 


11 


46.41 


11 


46.60 


11 


46.75 


11 


47.00 


11 


47.20 


11 


47.42 


11 


47.71 


11 


48.00 


11 


48.30 


11 


48.60 


11 


48.99 


11 


49.31 


11 


49.74 


11 


50.12 


11 


50.60 


11 


51.03 


11 


51.54 


11 


52.10 


11 


52.61 


11 


53.23 


11 


53.85 


11 


54.52 


11 


55.22 


11 


55.99 


11 


56.73 


11 


57.60 


11 


58.48 
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86.70 


5. 19800526E 


07 


6 .00219496E 


11 


59.42 


131 


26. 28 


5. 261092 89E 


07 


6.04892296E 


11 


60.42 


132 


-35.20 


5. 32460162E 


07 


6 . 09565096E 


11 


61.48 



RESERVOIR DEPLETED 



VITA 



Ossian R. Butterfield was born at Leominster, Massa- 
chusetts, on May 12, 1927, the son of Ossian R. and Helen P. 

Butterfield. He attended public schools in Leominster, 
graduating from Leominster High School in 1944. He entered 
the U. S. Navy in June, 1944, was enrolled in the Navy V-12 
College Training Program and attended Williams College in 
W i 1 1 i ams t o wn , Massachusetts. In November, 1945, he was 
transferred to the U.S. Naval Reserve Officers Training 
Course Unit at Brown University. In 1947, he was graduated 
from Brown with the degree of Bachelor of Science in Engi- 
neering. Upon graduation, he was commissioned an Ensign in 
the U.S. Navy Civil Engineering Corps. Duty assignments have 
included engineering duties at the Naval Proving Grounds, 
Dahlgren, Virginia; the Naval Air Station, Dallas, Texas; the 
Naval Stations in Argentia, Newfoundland, and in Bermuda, The 
West Indies; the Naval Air Station, Patuxent River, Maryland; 
the Naval Radio Station, Kamiseya, Japan; and the Naval Mine 
Warfare School, Yorktown, Virginia. In 1955 and 1956 he 
attended Rensselaer Polytechnic Institute in Troy, New York, 
where he received the Bachelor of Civil Engineering degree in 
June 1956. He has served as Executive Officer and as Command- 
ing Officer of U.S. Naval Mobile Construction Battalion Six 
deployed to Rota, Spain, and other bases in the Atlantic Ocean 



Area. His most recent duty assignment was in the Bureau of 
Yards and Docks, Navy Department, Washington, D.C., where he 
was Director of the Construction Division responsible for the 
management of projects engineered and constructed by the U.S. 
Navy for other governmental agencies. 

He was selected for postgraduate education in petroleum 
engineering and entered the Graduate School of The University 
of Texas in September, 1963. 

Butterfield was married to Ann Churchill of King George, 
Virginia, on April 11, 1950. They have two children, Marcia, 

age 13, and Brian, age 12. 

Permanent address: 6600 Willamette Drive 

Austin, Texas 
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